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Abstract  of  Dissertation  Presented  to  the  Graduate  Council 
in  Partial  Fulfillment  of  the  Requirements  for  the  Degree  of 

Doctor  of  Philosophy 

STUDIES  OF  ARTERIAL  BRANCHING  IN 
MODELS  USING  FLOW  BIREFRINGENCE 

By 

William  Joseph  Crowe,  Jr. 

March,  1969 

Chairman:   Dr.  R.  W.  Fahien 

Co-chairman:   Dr.  L.  J.  Krovetz 

Major  Department:   Chemical  Engineering 

The  optical  technique  of  flow  birefringence  (double  refraction) 
has  been  used  to  investigate  steady  and  pulsatile  flows  through  drilled 
Plexiglas  branches  simulating  those  of  the  arterial  system.   This 
technique  was  also  used  to  trace  flow  disturbances  created  by  inserting 
needles  and  cardiac  catheters  radially  and  axially  into  a  laminar  flow 
through  a  circular  tube. 

The  birefringent  liquid  used  was  a  1.4%  aqueous  suspension  of  the 
aniline  dye  Milling  Yellow,  and  birefringence  patterns  were  photographed 
on  16  mm  motion-picture  film.  The  experiments  covered  a  Reynolds-number 
range  from  0  to  250. 

A  scheme  has  been  proposed  for  calibrating  the  isochromatic-band 
distribution  obtained  for  steady  laminar  flow  through  a  circular  tube 
against  the  average  shear  stress  along  the  light  path  through  the  fluid. 
The  optical  effects  arising  from  curved  boundaries  and  differences  in 
refractive  index  have  been  included  in  this  calibration. 
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Analysis  of  birefringence  photographs  has  revealed  that  regions 
of  flow  stasis  and  boundary-layer  separation  occur  along  the  outer  walls 
of  branches  for  upstream  Reynolds  numbers  between  20  and  50,  depending 
on  the  angle  of  branch.   In  addition,  it  has  been  shox^ni  that  regions  of 
locally  increased  shear,  where  stresses  exceed  the  upstream  wall  stress, 
occur  along  the  inner  walls  of  branches  if   1.   the  branching  angle  is 
greater  than  30°,   2.   the  upstream  Reynolds  number  is  greater  than  150 
and   3.   the  downstream  arms  of  the  branch  are  of  approximately  the  same 
diameter.   Birefringence  films  for  pulsatile  flows  have  shown  that  these 
regions  do  not  move  with  the  incoming  pulse  but  remain  centered  approxi- 
mately one  diameter  downstream  from  the  point  of  bifurcation. 

It  has  been  suggested  that  locally  Inrreasprl  shp^rino  sfT-ess  nn 
the  arterial  intima  at  sites  of  branching  represents  that  physical 
mechanism  for  initiating  wall  trauma  and  atherosclerotic  disease  which 
is  most  closely  supported  by  pathologic  data. 

Photographic  evidence  has  been  presented  for  the  occurrence  of 
spiral  secondary  flows  in  branches,  and  the  importance  of  the  resultant 
mixing  effect  has  been  discussed  with  respect  to  the  Stewart-Hamilton 
technique  for  measuring  cardiac  output. 

The  propagation  of  flow  disturbances  has  been  characterized  by 
defining  branch  entrance  lengths  in  terms  of  distorted  birefringence 
patterns.   Entrance-length  data  and  birefringence  photographs  have  shown 
that  as  upstream  Reynolds  number  increases,  flow  disturbances  in  branches 
are  generated  and  propagated  by  the  following  three  processes:   1.   the 
readjustment  of  a  laminar  velocity  profile,   2.   the  onset  of  boundary- 
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layer  separation  and  3.   the  onset  of  secondary  flow. 

Both  Milling-Yellow  dye  and  water  have  been  used  to  study  the 
flow-distribution  characteristics  of  branches  for  steady  flows  in  the 
Reynolds-number  range  from  0  to  2100.   These  experimental  results  have 
shown  that  upstream  fluid  momenti'.m  per  unit  volume  is  the  principal 
dynamic  variable  affecting  branching  flow  distributions  and  that  area 
ratio  is  a  much  more  significant  geometric  parameter  than  is  branching 
angle. 

It  has  also  been  shown  that  wall  attachment  occurs  in  branches 
to  the  extent  that  side-arm  flow  will  exceed  that  in  the  main  downstream 
stem  providing  that  1.   the  branching  angle  is  less  than  60°,   2.   the 
upstream  Reynolds  number  is  below  a  certain  value  which  is  dependent  on 
the  fluid  and  3.   the  diameter  of  the  side  arm  is  approximately  equal 
to  that  of  the  main  downstream  stem. 
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SECTION  I 
INTRODUCTION 

During  the  past  twenty  years,  engineering  research  has  undergone 
a  transition  characterized  by  continual  erosion  of  once  well-defined 
boundaries  separating  individual  disciplines.   One  need  look  no  farther 
than  the  modern  electronics  industry,  with  its  increasing  dependence  on 
such  fields  as  metallurgy,  physical  chendstry  and  chemical  engineering, 
to  appreciate  this  fact.  What  might  be  called  the  interdisciplinary 
revolution   has  grown  both  in  scope  and  intensity,  affecting  not  only 
engineering  disciplines  but  also  creating  new  interfaces  with  the  life 
and  social  sciences. 

Undoubtedly,  several  factors  have  contributed  to  this  process 
but  perhaps  two  of  the  most  significant  have  been  first,  the  development 
of  new  materials  and  devices  having  broad  application  (e.g.,  the  high- 
speed electric  computer  and  the  laser)  and  second,  the  refinement  of 
existing  measuring  instruments  and  techniques,  coupled  with  the  ever- 
present  desire  of  researchers  to  extend  known  theoretical  and  experimen- 
tal results  to  new  areas  of  research. 

The  field  of  biomedical  engineering  has  evolved  as  a  direct 
consequence  of  this  process  and  is  now  enjoying  rapid  growth  both  in 
this  country  and  internationally.   Significantly,  this  growth  has  occurred 
at  both  academic  and  industrial  levels  and  is  reflected  by  new  curricula 
and  publications  in  the  field. 


Among  those  topics  receiving  attention  in  the  bioinedical- 
engineering  literature,  none  has  received  more  than  the  cardiovascular 
system.   This  is  not  surprising  when  one  considers  that  diseases  of 
the  circulation  are  among  the  most  common  yet-unsolved  problems  in 
medicine  and  that  the  circulatory  system,  featuring  unsteady-state  heat, 
mass  and  momentum  transport  within  materials  possessing  non-Newtonian 
rheology  and  functioning  under  a  closed-loop  automatic  control  system, 
involves  many  of  the  most  interesting  and  challenging  engineering  prob- 
lems of  the  day. 

The  great  majority  of  biomedical-engineering  studies  dealing  with 
the  cardiovascular  system  fall  into  one  or  both  of  the  following  catego- 
ries:  those  presenting  experimental  data  obtained  from  in  vivo^    in  vitro 
or  model  experiments,  and  those  proposing  r^tiiematical  models  and  analogs 
to  incorporate  existing  experimental  data.   The  studies  of  McDonald, 
Taylor,  Hardung,  Womersley,  Attinger,  Bergel,  Noordergraaf  and  Gessner 
are  a  few  examples  of  these. 

However,  with  some  notable  exceptions, ^^'^'^'^^'^^''^^''^^ '"^^'"^^'^^'^"^ 
few  studies  have  attempted  to  correlate  either  experimental  data  or 
theoretical  predictions  with  abnormal  clinical  and  pathological  findings. 
That  is  to  say,  little  attention  has  been  given  to  providing  physical 
bases  or  mechanisms  for  the  formation  of  observed  pathological  phenomena, 
as  opposed  to  providing  models  simply  to  explain  observed  experimental 
data.   The  recent  work  of  Fry^%emonst rating  that  excessive  shearing 
stress  causes  wall-tissue  erosion  is  a  particularly  significant  step  in 


this  direction. 

The  present  work  represents  an  attempt  to  further  this  approach 
by  focusing  attention  on  experiments  with  models  simulating  arterial 
branches  and  interpreting  experimental  results  in  the  light  of  known 
findings  from  the  medical  literature.   The  approach  here  is  then  in 
contrast  with  an  engineering-science  approach  to  the  problem,  wherein 
one  would  be  concerned  with  solving  the  equations  of  change  for  various 
types  of  fluids  and  wall  conditions. 

A.   Problem  Statement 

The  problem  may  be  posed  generally  in  two  parts:   First,  in  the 
flow  range  common  to  the  arterial  system  in  man,  what  physical  phenomena 
may  be  attributed  to  the  presence  of  a  branch  in  the  flowing  stream? 
Second,  what  are  the  implications  of  these  phenomena  with  respect  to  wall 
trauma  and  atherosclerotic  disease,  known  to  have  a  predilection  for 
formation  at  such  sites  in  the  arterial  system?  More  specifically,  the 
present  study  seeks  to  investigate  three  major  areas  of  branching  dynamics 
under  conditions  of  both  steady  and  pulsatile  flow  by  using  the  optical 
technique  of  flow  birefrigence.   These  three  areas  are:   1.   the  distri- 
bution of  wall  shearing  stresses  in  branches  insofar  as  determining 
whether  local  concentrations  might  occur  which  could  damage  the  vessel 
wall,   2.   branching -flow  distributions  as  functions  of  branching  angle 
and  upstream  Reynolds  number  and   3.   convective-mixing  effects  at 
branches  and  the  propagation  of  disturbances  downstream  as  these  relate 
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to  suspended  particles  in  the  blood  stream  and  to  the  Stewart-Hamilton 

19  57 
indicator-dilution  technique   '    for  determining  cardiac  output. 

With  the  exception  of  one  wye  branch,  this  study  has  been 
restricted  to  branches  of  the  sLngle- side-arm  type  fabricated  from 
drilled  or  milled  Plexiglas  block.   Though  model  branches  of  this  type 
cannot  imitate  branches  of  the  arterial  tree  in  every  detail,  particu- 
larly their  distensibility  and  wall  characteristics,  at  least  some 
indications  of  possible  flow  behavior  in  the  arterial  branches  can  be 
obtained  in  this  way;  and  for  those  branches  rigidly  bound  by  surrounding 
tissue  and  muscle,  such  a  model  represents  at  least  a  good  first  approxi- 
mation. 

The  approach  here  is  primarily  experimental  and  considers  only 
mechanical  interactions  between  the  flowing  stream  and  the  branch. 
Though  biochemical  processes  do  serve  a  vital  function  in  the  cardiovas- 
cular system,  there  is  no  reason  to  suppose  a  priori   that  biochemical 
activity  at  a  bifurcation  would  be  any  different  from  that  occurring  in 
unbranched  arterial  segments.   On  the  other  hand,  there  is  every  reason 
to  suppose  that  flow  behavior  at  two  such  sites  would  be  different,  per- 
haps even  to  the  extent  of  either  altering  normal  biochemical  function 
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resulting  in  the  development  of  arterial  disease,  as  Rodbard's  work 

suggests,  or  simply  causing  trauma  to  the  vessel  wall,  such  as  that 

shown  in  Fry's  work.    This  point  is  further  emphasized  by  the  studies 

22 

of  Texon,  Tmparato,  Helpern  and  Lord,  who  have  shown  experimentally 
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*   that  atherosclerosis  can  be  produced  artificially  in  dogs  by 

altering  blood-vessel  geometry. 

There  are  four  important  reasons  why  flow  birefringence  was 
chosen  as  the  flow-visualization  technique  used  in  the  present  study: 
1.   Since  flow  birefringence  is  an  optical  technique,  there  is  no  need 
to  introduce  probes  or  other  devices  into  the  flow  field  thereby  causing 
external  disturbances.   2.   It  is  possible  to  study  flow  behavior  through- 
out the  entire  branch,  not  just  along  a  single  streamline  as  with  the 
technirve  of  dye  injection.   3.   Prior  two-dimensional  studies-^^'^^'^^'''*! 
have  indicated  that  flow-birefringence  patterns  can  be  correlated  with  the 
shearing-stress  distribution,  so  the  possibility  of  obtaining  shear-stress 
information  in  branches  directly  seemed  fruitful.   4.   It  was  felt  that 
flow  birefringence  would  provide  a  convenient  way  of  tracing  flow  distur- 
bances and  entrance  effects  downstream. 

As  early  as  1954,  Wayland  had  suggested  the  possibility  of  using 
flow  birefringence  to  study  flow  disturbances  in  branches,^-*-  though  this 
was  unknown  to  the  author  until  recently.   Fortunately,  the  present 
study  has  incorporated  most  of  Wayland 's  suggestions. 

B.   Prior  Branching  Studies 
Interest  in  the  problem  of  branching  flows  is  not  new.   In  fact, 
problems  of  flow  distribution  have  been  of  engineering  concern  since 
ancient  times,  and,  today,  relationships  between  vessel  branching  and  flow 


behavior  are  receiving  considerable  attention  in  connection  with  the 
design  of  fluid  amplifiers  and  fluidic  digital-logic  devices. 

On  the  other  hand,  physiologic  interest  in  branching  has  arisen 
only  within  the  last  fifty  years  and  has  been  sporadic.   Recent 
studies,   '   demonstrating  that  atherosclerosis  shows  a  predilection  for 
formation  at  sites  of  branching  and  curvature  in  the  arterial  system, 
have  increased  the  interest  in  branching.   There  now  seems  little  doubt 
that  problems  associated  with  vessel  branching  on  a  small  scale  consti- 
tute significant  research  areas  both  in  engineering  and  medicine. 

Surprisingly  little  work,  either  theoretical  or  experimental, 
has  been  published  on  the  branching  problem.   Poisson's  work,^^  concern- 
ing reflected  and  refracted  waves  at  a  branch  site  was  among  the  earliest 
theoretical  studies  of  branching.   Poisson's  results  were  later  repeated 
by  Lord  Rayleigh  ^°  in  his  discussion  of  the  propagation  of  sound  through 
branched  tubes. 

The  earliest  physiologic  interest  in  vessel  branching  appears  to 
have  been  that  of  Moens^^  and  Frank, ^  who  were  concerned  with  the 
propagation  of  pulsewaves  through  the  arterial  tree.   Though  Frank's 
study  has  often  been  cited,  it  apparently  has  never  been  criticized, 
which  is  surprising  in  view  of  the  fact  that  it  contains  numerous  mathe- 
matical errors.   Frank  did  discuss  the  problem  of  a  pulsewave  encounter- 
ing a  bifurcacion  but  only  to  the  extent  of  repeating  Poisson's  formula- 
tions (incorrectly)  as  they  appeared  in  Rayleigh 's  work. 


More  recent  studies  of  branching  have  emphasized  'joth 
theoretical  and  experimental  approaches.   Womersley's  study^^  represent3 
the  most  comprehensive  theoretical  attack  on  the  problem  of  mathemati- 
cally modeling  pressure  and  flow  relationships  in  the  arterial  tree. 
His  proposed  model  was  based  on  the  linearized  Navier-Stokes  equations 
and  the  equations  for  an  ideally  elastic  wall.   Womersley  also  considered 
wave  reflections  at  branches  and  investigated  the  frequency  dependence 
of  wave-reflection  coefficients  and  input  impedence. 

Acrivos,  Eabcock  and  Pigford   Cvmsidered  the  problem  of  flow 
into  an  idealized  manifold  consisting  of  uniformly  spaced  side  ports  of 
equal  area.   They  proposed  a  simplified  matliematical  model  based  on 
friction  factors  and  pressure -recovery  coefficients  and  were  able  to 
show  good  agreement  with  their  pressure  and  flow  data  for  both  blowing 
and  sucking  manifolds. 

Knox   considered  the  more  general  problem  of  measuring  pressure 
drops  and  flow  distributions  for  branches  of  varying  branch  angle.   He 
considered  branches  of  the  single-side-arm  type  and  determined  pressure 
recoveries  at  bifurcations.   He  defined  a  recovery  length  for  branches 
as  the  distance  dotmstream  from  the  branch  necessary  for  the  pressure 
gradient  to  regain  linearity,  and  his  data  showed  that  recovery  lengths 
for  the  side  arms  were  significantly  greater  than  those  for  the  main 
stems. 

The  technique  of  dye  injection  was  used  by  Gopher  and  Dick,-*--^ 
Hahn,  Donald  and  Grier^^  and  by  Barnett  and  Cochrane^  to  study  particle 
and  flow  distribution  in  glass  models  simulating  venous  flow.   Distinct 


flow  channeling  was  noted  in  these  experiTients . 

20  21 
Helps  and  McDonald   '   investigated  venous-flow  patterns  both 

in  vivo   and  in  a  glass  model  and  found  that  at  Reynolds  numbers  less 
than  1000  injected  dye  formed  a  parabolic  profile  in  the  tributary, 
but  at  the  junction  circulating  movement  developed.   Two  sets  of  secon- 
dary flow  were  set  up  in  each  pipe  half,  and  this  persisted  "for  some 
distance  down  the  parent  trunk  until  a  new  parabola  was  formed  across 

the  single  trunk". 

33 

McDonald   used  the  dye- injection  method  to  show  that  marked 

disturbances  were  present  at  the  aortic  bifurcation  in  a  rabbit  at 
Reynolds  numbers  less  than  1000.   He  reported  that  dye  could  be  seen  to 
impinge  on  the  wall,  and  eddies  were  set  up  near  the  orifices  of  the 

two  branches. 

25  27  28 
Krovetz   '   '   conducted  similar  studies  using  glass  models  and 

35 
Silastic  (silicone  rubber)  molds  of  dogs'  arteries   to  simulate  arterial 

branching.   He  observed  and  photographed  secondary  flows  which  dis- 
tributed dye  particles  across  the  cross  section  as  if  the  flow  were 
turbulent,  even  though  the  upstream  flow  was  well  within  the  laminar 
range.   He  found  that  dye  particles  were  distributed  unevenly  between 
the  downstream  branches  if  the  incoming  dye  front  was  laminar  but 
were  distributed  evenly  if  the  incoming  front  was  turbulent.   He  also 
found  that  the  critical  upstream  Reynolds  num.ber  for  the  onset  of  secon- 
dary flow  in  branches  ranged  from  58  to  89%  of  the  critical  Reynolds 
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number  for  the  onset  of  true  turbulence  in  a  straight  tube. 

Another  study  by  Stehbens,  ^  reporting  a  series  of  experiments 
on  turbulence  in  glass  models  simulating  arterial  bifurcations  and  the 
curvature  of  the  carotid  artery,  found  critical  Reynolds  numbers  rang- 
ing from  306  to  1473  "calculated  for  the  lowest  rate  of  flow  at  which 

turbulence  could  be  induced  in  the  models". 

27 

Typical  experimental  results  obtained  by  Krovetz   are  shown 

in  Figures  1  and  2.   Figure  1  shows  a  laminar  dye  front  encountering 
a  branch  in  a  glass  tube,  and  Figure  2  is  the  equivalent  experiment 
for  flow  through  a  plastic  mold  of  a  dog's  artery.   In  both  cases, 
pronounced  secondary  flow  occurs  immediately  distal  to  the  bifurcation, 
and  significant  mixing  occurs  in  this  area.   In  the  glass  branch 
(Figure  1),  this  mixing  proceeds  downstream  approximately  five  diameters 
before  any  dye  enters  the  side  arm;  and  in  the  mold  (Figure  2) ,  no  dye 
at  all  enters  the  upper  downstream  branch,  indicating  a  definite 
channeling  effect  similar  to  that  observed  by  Barnett  and  Cochrane.-* 
These  figures  also  show  that  most  of  the  mixing  activity  occurs  near  the 
inner  walls  of  the  branches  (the  medial-crotch  region)  and  that  little, 
if  any,  mixing  occurs  near  the  outer  walls  in  the  immediate  region  of  the 

branches. 

-  14 

Other  studies  by  Attinger-^  and  by  Fox  and  Hugh   have  confirmed 

that  the  regions  along  the  outer  walls  in  branches  represent  areas  of 

boundary-layer  separation  and  local  stasis.   Attinger  used  suspended 

dust  particles  to  trace  flow  disturbances  in  branch  models  for  Reynolds 


10 


REYNOLDS       NUMBER       1440 


/ 


I  cm 


Figure  1.   Dye-streamline  experiment  (Evans  Blue  in  water)  by 
Krovetz  for  steady  flow  through  a  glass  branch.   Pronounced  secon- 
dary flow  and  mixing  occur  at  the  branch  even  though  upstream  flow 
is  laminar . 
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numbers  up  to  10,000.   In  addition  to  observing  helical  flow,  he  also 
observed  cavitation  (which  he  attributed  to  air  dissolved  in  the  ' 
stream)  at  the  highest  Reynolds  numbers  in  the  regions  of  boundary-layer 
separation  along  the  outer  walls.   Fox  and  Hugh  considered  open-channel 
flows  and  used  aluminum  powder  sprinkled  on  the  flowing  surface  to 
trace  disturbances  in  branched  configurations.   They  postulated  that 
the  boundary-layer  separation  and  stasis  along  the  outer  branch  walls 
are  contributing  factors  for  the  deposition  of  atheromatous  placque 
there. 


SECTION  II 


THEORETICAL  CONSIDERATIONS 


In  this  section,  we  present  the  basic  fluid-dynamic  and  optical 
theory  underlying  the  birefringence  technique  used  and  derive  the 
macroscopic  balances  of  mass  and  momentum  for  a  "generalized"  branching- 
flow  system.   Although  the  intent  here  is  primarily  to  develop  those 
mathematical  relations  needed  later  for  analyzing  the  flow  photographs, 
a  brief  discussion  of  flow  birefringence,  together  with  a  review  of 
certain  previous  studies  employing  this  technique,  has  been  included. 

A.   Fundamental  Equations  -•■-, 

We  start  by  considering  the  equations  of  continuity  and  motion 
in  their  Lagrangian  form  : 


Dg 

Dt  ~  -pV'v  continuity     (1) 


DV 
•^  Dt  =  -^P  -'  •  i  +  P£.  motion         (2) 


Here,  p  is  the  fluid  density,  v  the  fluid  velocity,  P  the 
pressure,  t^  the  shear-stress  tensor  and  _g^  the  vector  sum  of  body 
forces  per  unit  mass  acting  on  the  fluid.   D/Dt  is  the  substantial 
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derivative  defined  as  -r—    +     v 

0  C 


It  is  well  known  that  these  equations  apply  in  general  for  any 
continuous  medium  and  in  Lagran.jLan  form  apply  specifically  to  elements 
of  fluid  moving  with  the  flow. 

We  wish  to  consider  the  set  of  equations  resulting  from  (1)  and 
(2)  for  flow  through  a  uniform  cylindrical  tube  of  inside  radius  R. 
Anticipating  the  experimental  flow  conditions  to  be  encountered  later ^ 
we  use  the  cylindrical  coordinates  r,  G  and  z  and  impose  the  following 
assumptions  upon  Equations  (1)  and  (2) . 


^    .     0 

9t     - 


steady  flow 


=  constant 


3.     Vg   = 


4.  ^^ 

32 


0,  1 


3t 


rr 


rz 


=  0 


0 

0 
0 


^rz 

0 

^zz 


3t^ 


incompressible 
flow 


=  0^      angularly- 
symmetric  flow 


fully-developed 
f  lov7 


5.   g  =  0 


no  body  forces 


The  symmetry  of  the  stress  tensor  is  recognized  in  assumption  3,  as  is 
the  possible  existence  of  normal  stresses  (such  as  would  occur  in  a 
Stokesian  fluid). 
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From  the  continuity  equation  (1) ,  we  have  by  assumption  2 
that  V*  V  =  0.   That  is,  in  cylindrical  coordinates 


I     3(rvr)    1  av0      9 


V 


+-  — ^  +  — - —  =  0  (3)  -■■(■'^^ 


Since  the  second  term  of  Equation  (3)  is  zero  by  assumption  3 
and  the  third  term  zero  by  assumption  4,  the  continuity  equation 
becomes  simply 


linjl  =    0  (4) 

Integrating  Equation    (4)   directly,   we  obtain  '  '         '-'^ 


rvr  =  G  (5)  -■^ 

■■-V- 

where  the  function  of  integration  G  is  independent  of  r.  .   ■'.}:A 

If  Equation  (5)  is  solved  for  v  ,  then  it  becomes  apparent  that 
G  must  be  zero  on  the  physical  basis  that  the  velocity  component  Vv- 
cannot  be  infinite  along  the  axis  r  =  0.   We  therefore  conclude  that 
v^  itself  is  zero  and  that  the  above  assumptions  describe  a  parallel 
flow  with  v^  the  only  remaining  velocity  component. 

The  components  of  the  equation  of  motion  (2)  in  accordance  with 
assumptions  1  through  5  are 


T.-  ->, 
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M.  =  _  1  3(rTrr) 
3r    "  r   3r 


r-component     '  ■     (6) 


■S— .  =  0  0-component  (7) 


3P  _   1  9(rTj.z)  ^  fa\ 

=  — r— '-^—  z-component  (8) 

oZ        IT     oZ 


Equation  (7)  shows  that  the  pressure  P  is  independent  of  0  and 
therefore  dependent  only  on  r  and  z. 

If  Equation  (6)  is  differentiated  partially  with  respect  to  z 
and  the  order  of  differentiation  interchanged  on  the  right-hand  side, 
oi?e  has 


3_   iZ.  ^  ^1  3_  3(rTrr) 

3z   3r   "r  3r   8z 


(9) 


Since  the  right  side  of  (9)  is  zero  by  assumption  4,  it  follows  that 


1-   il  =  0  (10) 

3z   3r 


That  is,  the  radial  pressure  gradient  is  independent  of  axial  position. 
Because  the  cross  partial  derivatives  of  P  are  equal,  it  follows  that 


1-^=0  ■      .       (11) 

3r   37, 


That  is,  the  axial  pressure  gradient  is  independent   of  radial  position. 


'^'JjiB 
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According  to  Equation  (8) ,  the  axial  pressure  gradient  must  also  be 
Independent  of  both  0  and  z,  since  differentiation  of  that  equation 
with  respect  to  either  0  or  z  and  interchanging  the  order  of  differen- 
tiation on  the  right  side  renders  the  right  side  zero  by  virtue  of 

either  assumption  3  or  assumption  4  above.   Thus,-^  is  independent  of 

3z 
all  three  coordinates,  and  we  have  the  important  result  that  the  axial 

pressure  gradient  is  constant. 

Equation  (8)  may  now  be  integrated  directly  to  give 


r^rz  =  -  (.JT-^f-     +F  (12) 


where  the  function  F  is  independent  of  r. 

Solving  Equation  (12)  for  x   and  noting  that  t   cannot  be 

^2  rz 

infinite  at  r  =  0  leads  to  the  conclusion  that  F  must  be  zero,  and 
therefore 


'"  '  -  '!|)  r  (13, 


That  is,  the  shear  stress  ^^^  is  a  linear  function  of  radial  position. 

Equation  (13)  represents  the  starting  point  for  the  isochromatic 

shear-stress  calibration  to  be  discussed  in  Section  V.   There  we  shall 

find  it  convenient  to  drop  the  subscripts  on  r^^   and  to  consider 

Equation  (13)  in  terms  of  the  wall  stress  t  as  follows: 

w 


^  "  ^w-R  ■ (14) 
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where  t„  =  -  (  -^^  )  -y 


The  preceding  derivation  involved  integrating  the  equations  of 
continuity  and  motion  over  an  arbitrary  cross  section  of  a  circular 
tube.   We  now  consider  the  more  general  problem  of  integrating  these 
equations  over  an  arbitrary  volume  in  space  to  obtain  the  so-called 
macroscopic  balances  of  mass  and  momentum. 

This  problem  has  been  treated  by  Bird  for  a  system  consisting 
of  a  "generalized"  chemical  plant  having  a  single  inlet  flow  and  a 
single  outlet  flov?.   We  consider  the  equivalent  problem  for  a  "generalized" 
branching- flow  system  in  which  no  chemical  reactions  are  occurring  but 
which  may  have  any  number  of  inlet  and  cutlet  streams. 

The  arbitrary  region  of  interest  consists  of  that  volume  V 
bounded  laterally  by  the  walls  of  the  branch  and  bounded  at  the  inlets 
and  outlets  by  imaginary  control  surfaces  which  are  assumed  fixed  in 
space  and  time  and  oriented  perpendicular  to  the  direction  of  mean  flow 
through  the  inlet  or  outlet. 

It  is  convenient  to  start  from  the  equations  of  continuity  and 
motion  in  Eulerian  form: 


i^  =-  -  V  p  V  continuity         (15) 

3t  —      ■ 


— ti=  =  _V  •  pw  -Vp-7  •  i^  +  Pg  motion  (16) 

3t         —         -  -o. 
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These  equations  are  easily  derived  from  the  correspond in,<^  Lagrangian 
relations  (1)  and  (2),  but  when  written  in  the  above  form,  they  apply 
to  stationary  elements  of  volume. 

Integrating  the  continuity  equation  over  V  we  have  formally 


//I  If  '^V  =  --^^^  V  .  p  V  dV  (17) 

V  V 


By  an  extension  of  the  Leibnitz  rule  for  differentiating  an 

Q 

integral,   the  left  side  of  (17)  may  ho  written  as 

//I  |P  dV  =  f^  I//  P  dV-/J  pv^  .  n  dS  (18)   . 

V  V      s 

where  S  refers  to  the  surface  bounding  the  branch,  v  is  the  velocity 
of  this  surface  at  any  point  and  n  is  the  outwardly  directed  unit  vector 
normal  to  the  surface.   Note  that  Vg  is  zero  for  a  rigid  branch,  but 
would  not  be  so  in  the  cardiovascular  system. 

By  the  Gauss  divergence  theorem,  the  right-hand  integral  of  (17) 
becomes 


(jj   V-  p  v  dV=  //  p  v-  n  dS  (ig-) 

V  S      ~ 


where  the  surface  integral  is  again  taken  over  the  bounding  surface  S, 
and  v  is  the  velocity  of  the  fluid. 

Substituting  Equations  (18)  and  (19)  into  (17)  and  rearranging 
terms  gives 
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^  //^  P  dV  =  -//  p  rv  -Vg)  .  n  dS  (20) 

The  left  side  of  Equation  (20)  is  the  time  rate  of  change  of  the 
total  mass  M  within  the  branch.   The  integrand  on  the  right  side  vanishes 
on  all  solid  surfaces,  and  v^  vanishes  on  the  control  surfaces  at  all 
Inlets  and  outlets.   If  we  now  assume  that  the  density  p  does  not  vary 
across  the  inlet  and  outlet  surfaces,  and  note  that  the  normal  n  is 
opposite  in  sense  to  v  for  inlet  streams  but  in  the  same  sense  as  v 
for  the  outlet  streams  we  may  carry  out  the  above  integration  to  obtain 

dM   _  r-->   —  c      .— .   — 

dT  -Epi-i^i  -E^j-jS.  (21) 

IN         OUT 
where  the  first  summation  is  taken  over  all  inlet  streams  and  the  second 
over  all  outlet  streams.   The  subscripts  i  and  j  refer  to  the  ith  inlet 
stream  and  jth  outlet  stream  respectively.   S.  and  S.  are  the  cross- 
sectional  areas  of  these  streams,  and  the  average  velocity  v  is  by  definition 

J/^vdS 
v  = 


//dS 
S 


where  v  is  the  magnitude  of  v  at  the  cross  section  of  interest. 

If  we  now  define  the  difference  operator  A  as  meaning  the  sum  over 
all  inlet  streams  minus  the  sum  over  all  outlet  streams,  then  Equation 
(21)  may  be  conveniently  retvritten  as 


dM  _    ,/  -  ^^ 

d^  -  -A(P  V  S)  (22) 


•  fr^'ma 
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Equation  (22)  is  the  macroscopic  or  over-all  unsteady-state 

mass  balance  for  the  branching-flow  system.   Under  steady-flow 

conditions,  dM  =  o. 
Zt 

Integrating  the  equation  of  motion  we  have 

///  4^^^'  dV  =  -///(V  •  P  Z  v)dV  -///  (vP)dV  -///  (V.  T^)dV    (23) 
V  V  V  V 

+  //I  (pg.)dV 
V 

The  left  side  of  this  expression  may  be  transformed  as  before 

to  give 

r/rlMldV  =  ^|//(Pv)dV  -/|(Pv)v_^  •  n  dS  (24) 


The  first  three  integrals  on  the  right  side  of  (23)  may  be  transformed 
to  surface  integrals  by  modifications  of  the  divergence  theorem  as 
follows  : 


///  (V-  P  V  v)  dV  =  /f  (p  v)v  •  n  dS  (25) 

V  S 


jjj    (V  P)dV  =   1/  P  n  dS 


(26) 


///  (V-  X  )  dV  =  /J  i  •  n  dS  C27) 

V  s 
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Substituting  Equations  (24)  through  (27)  into  (23)  and 
rearranging  terms  we  have 

—  //I  p  V  dV  =  -//  p  V  (v-Vg)  •  n  dS  -//  P  n  dS  (28) 

V         s  s 

-/I  (X  •  n)dS  +  ///   p  ^  dV 
S  V 

The  left  side  of  Equation  (28)  is  the  time  rate  of  change  of 
the  total  momentum  P^  within  the  branch.   Again,  the  first  integrand  on 
the  right  side  vanishes  on  all  solid  surfaces,  and  v  vanishes  at  all 
inlets  and  outlets.   The  remaining  contribution  from  this  term  points 
in  the  direction  of  the  flow  both  at  the  inlets  and  outlets.   The  second 
term  on  the  right  contributes  to  the  influx  and  efflux  of  momentum  at 
the  inlets  and  outlets  respectively  and  to  the  drag  force  on  all  solid 
surfaces.   The  third  term  represents  the  viscous  contribution  to  drag  on 
solid  surfaces  and  gives  a  usually  negligible  contribution  to  momentum 
transport  at  flow  inlets  and  outlets.   Finally,  the  last  term  on  the 
right  represents  the  vector  sum  of  all  the  body  forces  F^odv  acting  on 
the  fluid. 

Using  the  previously  defined  A  notation,  we  may  write  the  macroscopic 
momentum  balance  for  the  branching-flow  system  as 


dP  "T 

-~r     =     -  A  (p  v'^  S  +  PS )  -•  F,    +  F.  ,  (29) 

dt  —      _    ^j.gg   —body  ^  ^ 


'  ;  :-'*-'^'s«;"*"'< 
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where  the  vector  S^  has  a  magnitude  equal  to  the  cross-sectional  area 
of  the  inlet  or  outlet  of  interest,  but  points  in  the  direction  of  the 
flow.   F<jra£  ^"'^  '^  ^^^   defined  as  follows: 


Idrag  =  //  P  S.  dS  +  1/  1^  •  n  dS 
^solid      ^solid 


//cV^ 


^   = 


JJdS 

s 


Again,  for  steady-flow  conditions  -=  =  0. 

dt     ■" 


B. 
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According  to  Jerrard,   "the  first  published  observations  on 

30 
flow  birefringence  were  made  in  1873,  when  Mach   observed  double 

refraction  in  extremely  viscous  substances  such  as  strong  metaphos- 

phoric  acid  and  Canada  balsam,  when  poured  into  a  beaker."  Soon  there- 

31 
after.  Maxwell   reported  that  he  had  experimented  with  birefringence 

as  early  as  1866  and  had  devised  a  concentric -cylinder  apparatus  for 

studying  the  phenomenon  in  Canada,  balsam.   Since  that  time,  flow 

birefringence  (also  called  streaming  birefringence  and  flow  double 

refraction)  has  been  observed  in  many  liquids,  and  numerous  theories 

of  flow  birefringence  have  been  proposed.   From  the  extensive  reviews 

23  9 

of  these  theories  presented  by  Jerrard   and  by  Cerf  and  Scheraga, 

it  is  evident  that  no  single  theory  has  been  successful  in  predicting 
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all  birefringence  phenomena  for  liquids. 

Though  most  of  the  early  studies  of  flow  birefringence  were 
directed  toward  using  this  as  a  technique  for  investigating  molecular 
size  and  shape,  within  the  last  forty  years,  many  investigators  have 
reported  using  biref inngence  techniques  to  study  the  properties  of 
flows  themselves.   Attinger's  recent  studies"^'^  of  eddy  formation  and 
flow  disturbances  both  in  distensible  tubes  and  in  rigid  tubes  of  circu- 
lar and  elliptical  cross  section,  are  typical  examples  of  the  qualitative 
use  of  flow  birefringence.   The  most  significant  quantitative  work  has 
been  that  by  Prados^-  and  Prados  and  Peebles, ^°  who  showed  that  the  bire- 
fringence technique  could  be  used  to  determine  shear-stress  distributions 
and  velocity  profiles  for  flows  through  various  types  of  flat  channels. 
They  used  aqueous  suspensions  of  Milling-Yellow  dye  and  obtained  good 
agreement  between  the  flow  rates  calculated  from  birefringence  measurements 
-and  the  experimentally  determined  values.   Essentially  no  quantitative 
work  has  been  attempted  for  three-dimensional  flow  systems,  and, apparently, 
the  first  attempt  to  use  flow  birefringence  for  studying  unsteady  or  pul- 
satile flow  was  the  1957  work  of  Thurston  and  Hargrove  cited  by  Prados  and 
Peebles.    It  is  probable  that  the  present  study  represents  the  first 
attempt  to  use  flow  birefringence  for  studying  pulsatile  flows  in  branches. 

As  is  the  case  for  photoelastic  solids,  flow  birefringence  in 
liquids  results  from  the  optical  anisotropy  induced  by  a  shearing  field. 
The  phenomenon  is  most  pronounced  in  solutions  or  suspensions  containing 
asymmetric  particles,  although  it  has  been  observed  at  high  shear  rates 
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in  pure  liquids  (and  Newtonian  fluids) ,  where  the  anisotropy 
presumably  arises  because  of  molecular  deformation  rather  than  molecu- 
lar reorientation. 

A  typical  flow-birefringence  experiment  requires  that  an  incident 
beam  of  plane-polarized  light  (from  the  polarizer)  be  passed  through 
the  flowing  liquid  in  a  direction  perpendicular  to  the  direction  of 
flow.   The  observer  then  "analyzes"  the  emerging  light  by  viewing  this 
light  through  a  second  polarizer  (the  analyzer)  which  is  oriented  so 
that  its  plane  of  polarization  is  crossed  90°  with  respect  to  that  of 
the  incident  light.   Tt  is  found  that  when  the  liquid  is  at  rest,  the 
incident  beam  of  polarized  light  emerges  unchanged  from  the  liquid  and 
is  extinguished  by  the  analyzer.   Thus,  the  entire  field  appears  dark 
to  the  observer  under  zero-flow  conditions.   However,  when  the  liquid 
is  set  in  motion,  the  observer  sees  a  series  of  interference  bands  dis- 
tributed through  the  flow  field.   These  bands  appear  dark  when  monochro- 
matic light  is  used  and  colored  when  white  light  is  used,  and  their 
presence  indicates  that  only  under  certain  conditions  does  the  incident 
beam  emerge  unchanged  from  the  flowing  liquid. 

Although  the  above  description  pertains  to  the  use  of  plane- 
polarized  light,  a  similar  set  of  results  is  obtained  when  circularly 
polarized  light  is  used.   An  incident  beam  of  circularly  polarized 
light  can  be  produced  by  inserting  a  quarter-wave  plate  between  the 
polarizer  and  the  liquid  such  that  the  optic  axis  of  the  quarter-wave 
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plate  makes  an  angle  of  45°  with  the  plane  of  polarization  of  the 
polarizer.   The  effect  of  this  quarter-wave  plate  is  to  decompose 
the  incident  plane-polarized  vibration  into  two  components  vibrating 
at  right  angles  to  one  another  and  90°  out  of  phase.   These  are  pre- 
cisely the  properties  of  circularly-  polarized  light.   \Jhen  circularly 
polarized  light  is  used,  it  is  also  necessary  to  insert  a  quarter-wave 
plate  between  the  fluid  and  the  analyzer.   This  plate  is  oriented  so 
that  its  optic  axis  is  crossed  90°  with  respect  to  that  of  the  first 
quarter-wave  plate.   If  a  light  beam  emerges  unchanged  from  the  fluid 
(that  is, circularly  polarized),  the  second  quarter-wave  plate  changes 
the  vibration  back  to  a  single  plane-polarized  beam  vibrating  perpen- 
dicular to  the  plane  of  polarization  of  the  analyzer.   The  bean  is  then 
extinguished.   Regardless  of  the  type  of  polarized  light  used,  it  is 
true  that  the  observer  will  see  an  interference  band  in  the  flow  field 
only  at  positions  where  the  incident  light  emerges  unchanged  from  the 
liquid. 

The  effects  of  passing  beams  of  plane  and  circularly  polarized 
light  through  doubly  refracting  materials  and  the  conditions  under  which 
such  beams  will  emerge  unchanged  have  been  treated  analytically  by 
Rosenberg.  ^     Basically,  if  a  beam  of  plane-polarized  light  enters  a 
birefringent  liquid,  this  beam  will  be  decomposed  into  two  component 
vibrations,  one  vibrating  parallel  to  the  optic  axis  of  the  fluid  and 
the  other  perpendicular  to  this  axis.   Initially,  these  two  vibrations 
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are  in  phase  and  may  be  represented  by 


X  =  a  cos  (277^)  (30) 


y  =  b  cos(27t1)  (31) 


where  x  is  the  displacement  of  the  component  vibration  parallel  to  the 
optic  axis  of  the  fluid,  y  the  displacement  of  the  component  vibration 
perpendicular  to  this  axis,  a  the  amplitude  of  the  x-vibration,  b  the 
amplitude  of  the  y-vibration,  and  T  is  the  period  of  vibration. 

Owing  to  the  optical  anisotropy  of  the  flowing  liquid,  these  two 
vibrations  will  in  general  propagate  through  the  liquid  at  different 
speeds  and  will  emerge  out  of  phase  by  an  amount  dependent  upon  the  path 
difference  (or  relative  retardation)  6  between  them.   The  emerging 
vibrations  may  therefore  be  expressed  as 


X  =  a  cos(2Tit)  (32) 

T 


y  =  b  cos  27r(^  +  |)  (33) 


where  X  is  the  wavelength  of  the  incident  beam. 

If  the  emerging  components  are  now  recombined  by  eliminating  t, 
the  resulting  equation  is 


x^  +  yi   _  2xy  cos  (2Tr^)       .  2  .   6^  .   . 

I?    b?         ab  =  """  ^2"A^  ^^^ 
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Since  Equation  (34)  describes  an  ellipse,  the  emerging  beam  will 
in  general  be  eiliptically  polarized.   However,  if  one  compares  Equa- 
tion (34)  with  Equations  (30)  and  (31),  it  is  apparent  that  under  the 
following  conditions  the  beam  will  emerge  unchanged 

6  =  NX,  N  =  0,  1,  2....  (35) 

a  =  0  or  b  =  0  (36) 

The  first  of  these  conditions  governs  the  appearance  of  the 
so-called  isochromatic  bands  (or  isochromatics) ,  while  the  second  con- 
dition specifies  the  appearance  of  the  so-called  isoclinics.   Thus,     ■■ 
when  plane-polarized  light  is  used,  dark  fringes  can  appear  at  the  analy- 
zer from  both  isochromatics  and  isoclinics. 

Equations  (35)  and  (36)  show  that  the  isochromatic  bands  are 
related  to  the  magnitude  of  the  path  difference  6,  while  the  isoclinics 
appear  only  when  the  incident  plane-polarized  vibration  is  parallel 
(b  =  0)  or  perpendicular  (a  =  0)  to  the  optic  axis  of  the  fluid.   It 
has  been  found  experimentally  that  the  path  difference  5  is  dependent 
upon  the  length  of  the  light  path  and  upon  the  average  shearing  stress 
along  this  path  while  the  orientation  of  the  optic  axis  is  dependent  upon 
the  direction  of  shearing  stress  along  the  light  path.   Thus,  it  is  clear 
that  the  isochromatics  are  related  to  both  the  magnitude  of  the  shearing 
stress  and  the  length  of  the  light  path,  while  the  isoclinics  are  related 
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to  the  direction  of  shearing  stress  along  the  light  path. 

The  analysis  is  similar  for  circularly  polarized  light. 
However,  since  circularly  polarized  light  consists  of  two  mutually 
perpendicular  vibrations  which  are  of  equal  amplitude  but  90°  out  of 
phase,  each  of  these  vibrations  will  be  decomposed  into  components 
vibrating  parallel  and  perpendicular  to  the  optic  axis  of  the  fluid. 
When  the  appropriate  contributions  from  the  incident  beam  are  added, 
the  components  propagating  through  the  fluid  may  be  represented  by 


X  =  a  cos  2tt^  (37) 


y  =  a  sin  21?^'  (38) 

T 


where  a  is  the  amplitude  of  the  incident  beam,  T  the  period  of  vibra^- 
tion  and  t'  =     t     +    —  with  e  being  the  angle  of  rotation  between  the 
abscissa  describing  the  incident  beam  and  the  optic  axis  of  the  fluid. 

Again,  these  components  will  undergo  a  relative  retardation 
resulting  in  a  corresponding  phase  change.   Upon  emerging  from  the  fluid 
the  components  can  be  represented  by 


X  =  a  cos  2Tpt  (39) 


y  =  a  sin  2^(l!-  +  i)  (40) 

1        A 
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Equations  (39)  and  (^0)  can  be  combined  by  eliminating  t'  to 
give 


x2  +  y2  -2xy  sin(2Ti-)  =  a^  cos  (2ttA)  (41) 

X  X 


Equation  (41)  also  describes  an  ellipse,  and, hence,  when 
circularly  polarized  light  is  used,  the  emerging  beam  will  in  general 
be  elliptically  polarized.   However,  comparison  of  Equation  (41)  with 
Equations  (37)   and  (38)  shows  that  the  circularly  polarized  beam  will 
emerge  unchanged  only  if  6  =  NX.   Therefore,  when  circularly  polarized 
light  is  used,  only  the  isochromatic  bands  can  appear  in  the  flow  field. 

It  has  been  shown  that  the  condition  5  =  NX  corresponds  to  con- 
structive interference  of  the  emerging  vibrations  for  both  plane  and 
circularly  polarized  light.   Since  the  two  vibrations  are  known  to 
travel  at  different  speeds  through  the  fluid,  constructive  interference 
also  implies  that  the  difference  between  the  number  of  waves  of  each 
vibration  in  the  fluid  is   an  integer  N.   Because  the  difference  in 

wave  speeds  of  the  two  vibrations  occurs  from  changes  in  wavelength  and 

52 
not  changes  in  frequency,   constructive  interference  requires  that 


(£/Xx  -a/Xy)      =  N  (42) 


where  Z   is  the  length  of  the  light  path,  X^  the  wavelength  associated  with 
the  vibration  parallel  to  the  optic  axis  and  X   the  wavelength  associated 
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with  the  vibration  perpendicular  to  the  optic  axis. 

Since  the  frequency  of  vibration  f  is  unchanged  in  the  fluid, 
the  wave  speeds  must  obey  the  following  equations 

Cx  =  f^x  (43) 

Cy    =    fAy  (44) 

The  wave  speed  of  any  light  ray  in  vacuo   is  given  by 

c  =  fX      ■  (45) 

where  c  is  the  speed  of  light,  f  the  frequency  of  vibration  and  \   the 
wavelength  of  the  vibration  in  vacuo. 

If  Equation  (45)  is  now  divided  by  Equations  (43)  and  (44)  the 
following  relations  are  obtained: 

A/Ax  =  c/c^  =  n^  (46) 


X/Xy  =  c/Cy  =  Uy  (47) 


where  n^^  and  ny  are  by  definition  the  indices  of  refraction  for  the  two 
vibrations. 

Equation  (42)  may  now  be  rewritten  in  terms  of  the  indices  of 
refraction  n^  and  n^  as  follows: 

(nx  -ny)  jl   =  NX  .  (48) 
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This  expression  gives  the  important  relation  between  the  so-called 
amount  of  birefringence    (n  -n^ )  and  the  isochromatic-band  order  N.  The 

X.  y 

quantities  9n^   and  ftiy  are  know  as  the  optical  path  lengths   for  vibra- 
tions X  and  y  respectively. 

According  to  Equation  (48)  the  phase  difference  in  radians  between 
the  two  vibrations  is  given  by 

A  =  2^(n^-ny)  j  ^^^^ 

A 

This  is  the  result  given  without  derivation  by  Jerrard.^^ 

The  quantitative  application  of  flow  birefringence  begins  with 
Equation  (48).   The  usual  assumption  employed  states  that  the  amount  of 
-birefringence  (n^  -ny)  is  a  linear  function  of  the  average  shear  stress 
T(avg)  along  the  light  path.   The  results  obtained  by  Prados^°  indicate 
that  this  is  a  good  assumption  for  suspensions  of  Milling-Yellow  dye. 
This  assumption  may  be  stated  mathematically  as 

(n^  -  ny)   =  MT(avg)  (50) 

where  the  coefficient  M  is  known  as  the  Maxwell  constant. 

If  the  assumption  given  by  Equation  (50)  is  substituted  into 
Equation  (48)  one  has 

NX  =  M£T(avg)  •  (51) 

Defining  M*  as  M/X,  and  rewriting  Equation  (51)  in  terms  of  the 
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Isochromatic-band  order  N  we  have  finally 

N  =  M*£T(avg)  (52) 

This  expression  represents  the  starting  point  for  the  analysis 
of  the  flow  photographs  to  be  presented  in  Section  V. 


H^'-r- 


SECTION  III 


EXPERIMENTAL  EQUIPMENT 


The  production  of  flow-birefringence  patterns  requires  a 
relatively  simple  experimental  setup.   Basically,  the  birefringent 
liquid  must  be  caused  to  flow  through  a  transparent  test  section 
located  between  two  crossed  polarizing  filters  (it  is  common  practice 
to  refer  to  the  incident-light  filter  as  the  polarizer  and  to  the  second 
filter  as  the  analyzer).   A  light  source,  which  may  or  may  not  be  mono- 
chromatic, must  be  oriented  such  that  the  incident  unpolarized  light 
passes  through  the  polarizer,  test  section  and  analyzer  respectively. 
The  observer,  looking  through  the  analyzer  in  the  direction  of  the  light 
source,  sees  optical  interference  patterns  wherever  shear  stresses  are 
present  and  a  dark  field  elsewhere. 

The  experimental  setup  used  in  the  present  study  was  designed  with 
simplicity  and  ease  of  operation  in  mind  and  is  pictured  schematically 
in  Figure  3.   For  purpose  of  discussion,  this  setup  may  be  conveniently 
subdivided  into  three  major  components:   A.   The  fluid-perfusion  system, 
B.   The  flow  test  sections  and  C.   The  optical  system.   Where  pertinent, 
the  important  pieces  of  peripheral  equipment  are  also  discussed. 

A.   The  Fluid-Perfusion  System 
Since  flow  birefringence  occurs  under  dynamic  conditions. 
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provision  must  be  made  for  either  recirculating  the  birefrigent  fluid 
through  a  test  section  or  collecting  it  after  a  single  pass.   It  was 
decided  that  for  branching-flow  studies,  and  in  particular  for  deter- 
mining branch  flow  ratios,  the  best  approach  was  to  provide  for  a 
single  pass  through  the  test  sections  with  fluid  collection  and  metering 
at  the  outlets.   Direct  unobstructed  flow  metering  using  graduated  cylin- 
ders was  thus  possible,  assuring  the  experimenter  that  observed  flow 
ratios  were  due  to  branch  configuration  rather  than  side  effects  of 
flow  metering.   The  single-pass  approach  provided  the  additional  advan- 
tage of  eliminating  the  need  for  a  constant-temperature  bath  to  remove 
viscous  heat  caused  by  continuous  recirculation.   Change  in  fluid  vis- 
cosity due  to  viscous  heating  was  a  difficulty  reported  by  Peebles, 

37 
Prados  and  Honeycutt  in  a  previous  study. 

General  requirements  were  that  the  fluid-perfusion  system  be 
-capable  of  producing  both  steady  and  pulsatile  flows  over  the  same  range 
of  flow:.rates,  that  the  fluid  be  maintained  at  constant  temperature  and 
isolated  from  outside  vibration  during  experimental  runs  and  that  the 
flow  system  be  leveled  so  as  to  eliminate  gravitational  effects  in  the 
flow  test  sections. 

Steady  flows  were  generated  using  the  flow  system  illustrated 
in  Figure  3.   This  system  consisted  of  a  glass  tank  approximately  30  cm. 
high  which  could  be  sealed  with  a  rubber  stopper  and  pressurized  using 
an  external  compressed-air  supply.   Though  a  larger  constant-head  tank 
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was  available  In  the  laboratory,  accurate  flow  control  over  a  much 
wider  range  of  flow  rates  was  possible  using  the  pressurized  system. 
A  pressure  gauge  permitted  tank  air  pressure  to  be  monitored  to  better 
than  1  mm  of  mercury  over  the  range  0  to  200  mm  of  mercury. 

The  fluid  tank  was  equipped  with  a  tapering  outlet  nozzle 
converging  to  an  outside  diameter  of  9  mm,  just  sufficient  to  match 
the  outside  diameter  of  the  glass  connecting  tube  joining  the  fluid 
tank  to  the  flow  test  section.   This  provided  for  a  smooth  transition 
flow  from  tank  to  connecting  tube.   The  connecting  tube  used  for  all 
experiments  was  a  3.00  mm  diameter  section  of  precision-bore  glass 
tubing,  96  cm  in  length.   This  length  was  sufficient  to  insure  that 
fully  developed  flow  entered  the  test  sections. 

A  tight-fitting  rubber  sleeve  over  the  nozzle  and  connecting 
tube  served  both  as  a  seal  against  leaks  and  a  structural  support,  thus 
assuring  a  secure  flush  connection  between  nozzle  and  connecting  tube. 
In  addition,  this  sleeve  provided  a  simple  means  for  adapting  the  flow 
system  from  steady- tD  pulsatile-flow  operation.   To  do  this,  it  was 
simply  necessary  to  slide  the  rubber  sleeve  off  the  connecting  tube 
while  pinching  it  to  arrest  any  flow  due  to  fluid  in  the  tank.   The 
glass  tube  could  then  be  connected  directly  to  the  pulsatile-flow  source 
using  a  section  of  rubber  tubing. 

Pulsatile  flows  simulating  heart  action  were  generated  using  a 
positive-displacement  rotary  bellows  pump  (Model  1000,  manufactured  by 
Research  Appliance  Company,  Pittsburgh,  Pennsylvania).   This  pump  ' 
employed  a  rotating  eccentric  cam  to  displace  fluid  through  a  ball  valve. 
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The  resulting  flow  wave  was  measured  independently  using  a  square-wave 
electromagnetic  flowmeter  (Model  301,  manufactured  by  Carolina  Medical 
Electronics  Company,  Winston  Salem,  North  Carolina)  and  is  shown  in 
Figure  4.  •  • 

Whether  operating  under  steady-or  pulsatile-flow  conditions,  the 
entire  fluid -perfusion  system  was  mounted  on  sturdy  laboratory  tables 
and  maintained  in  an  air-conditioned  laboratory  at  a  constant  ambient 
temperature  of  25°C.   External  vibration  was  small  but  could  be  mini- 
mized by  using  felt  pads  for  supports  and  by  turning  off  the  air-condi- 
tioning system  during  experimental  runs. 

B.   The  Flow  Test  Sections 

There  are  several  desirable  properties  for  flow  test  sections 
to  be  used  in  birefrigence  work.   Among  these  are  that  they  should  be 
as  nearly  transparent  as  possible  and  of  reasonable  cost.   They  should 
also  be  easily  fabricated  from  a  material  in  which  residual  stresses 
can  be  removed.   If  the  last  requirement  is  not  met,  the  resulting  defor- 
mation of  crystal  planes  inthe  material  may  render  it  optically  active 
with  respect  to  polarized  light,  hence  interfering  with  experimental 
results. 

While  there  may  be  no  ideal  material  in  every  respect,  previous 

40  41 
studies   '   have  reported  success  using  Plexiglas  test  sections.   This 

material  can  be  fabricated  with  relative  ease  and  is  reasonably  priced. 
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although  the  cost  of  labor  for  fabrication  may  be  considerable. 
Plexiglas  is  optically  active  in  polarized  light  when  residual  stress 
is  present,  but  this  can  be  removed  by  proper  temperature  curing.   Such 
test  sections  can  be  polished  both  inside  and  out  to  a  high  degree  of 
transparency,  rendering  them  ideal  from  the  optical  point  of  view. 
The  flow  test  sections  used  in  the  present  study  have  been 
listed  in  Table  1.   With  the  exception  of  sections  1,  2  and  7,  all  test 
sections  were  constructed  of  drilled  Plexiglas  block  and  were  fabricated 
in  various  branched  configurations  in  an  attempt  to  simulate  typical 
branches  found  in  the  cardiovascular  system. 

Fabrication  cost  imposed  a  hard  limit  on  the  number  and  variety 
of  branches  which  could  be  investigated. 

Figure  5  is  a  photograph  of  the  branched  test  sections  listed  in 
Table  1.  Over-all  size  of  these  branches  was  limited  by  the  lengths  of 
the  drills  available.  Test  section  7  was  milled  to  a  3  mm  wide  by  2  mm 
deep  rectangular  cross  section  so  that  an  experimental  comparison  could 
be  made  regarding  the  effects  of  internal  geometry  upon  optical  and  flow- 
ratio  behavior. 

Great  care  was  required  in  drilling  the  Plexiglas  block,  as 
eccentric  motion  of  the  drill  caused  internal  burring  which  could  not 
be  removed  by  later  polishing.   Care  was  also  necessary  to  insure  that 
branches  met  at  the  proper  angles  and  lay  in  a  plane  paralleling  that  of 
the  broad  face  of  the  Plexiglas  block. 

The  upstream  faces  of  all  blocks  were  drilled  to  a  depth  of 


o 
o 

o 


0) 

rH 

(0 

CO 

iH 

0] 

CO 

4J    -H 

« 

<0 

cu  u 

iH 

.H 

0)  -o 

O 

o 

o 

(U 

0) 

(U    O 

K 

V4 

Q 

o 

CD    Q) 

1 

1 

0)    k 

c 

c 

o         cfj 

o 

o 

C  r^    CO 

•H 

•H 

n       iH 

cn 

CO 

)-i    Jj    oo 

•H 

•H 

^     <U    -H 

o 

U 

-o    >< 

(U 

01 

rH     B     OJ 

v< 

>-i 

>-H    3   rH 

a< 

a. 

<   Z    CM 

41 


60 

c 
o 


1 

u 

0) 

(0 

a 

a 

o 

)-i 

o 

y-N 

t3 

t-l 

<U 

to 

.H 

>-t 

.-1 

3 

•H 

60 

^ 

C 

to 

4J 

C 

o 

o 

OJ 

•H 

(X3 

4-> 

0) 
60 
C 
O 
•J 


0) 
60 

c 
o 


to 

H 


o 

U 
W 
CO 

H 

W 

O 


CO 

(U 

u 

<u 

e 

CO 
•H 

O 


g 

O 

• 

•a 

•H 

o 

<u 

c 
o 

o 
o 

• 

O 
O 

• 

o 
o 

rH 

o 

o 
o 

o 
o 

to 

z 

z 

fO 

CN 

r-i 

CO 

CO 

fO 

CO 

• 

CO 

• 

CO 

g 

o 

• 

o 

4-1 
CO 

00 

CNJ 

c 

•H 
CO 

o 
o 

• 

o 

• 

o 
o 

O 
O 

o 
o 

1^ 

iH 

o 

o 

o 

o 
o 

rH 

E 

CO 

CO 

ro 

ro 

CO 

CO 

• 

• 
CO 

• 

CO 

• 

CO 

• 

CO 

e 

o 

• 

o 

to 

0) 

fO 

eg 

_ 

>-< 

X 

o 

4-1 
CO 

O 

o 

o 

o 
o 

O 
O 

o 

o 

o 
o 

o 
o 

3 

en 

ro 

ro 

CO 

CO 

CO 

CO 

• 
CO 

CO 

• 
CO 

• 

CO 

6( 
C 


s 


C 
O 


01 

c 
o 


o 

ro 


0 

0 

o 

o 

O 

o 

O 

o 

CO 

CO 

CO 

CO 

o 

o 

o 

o 

O 

O 

O 

m 

vO 

SO 

C7\ 

?l 

« 

CO 

o 


j= 

o 

c 

to 

h 

w 

B 

M 

0) 

CO 

I 

3 

(1) 

H 

,£ 

T) 

3 

•H 

S^ 

H 

to 

to 

1 

iH 

(U 

0) 

3 

)^ 

iH 

U 

to 

60 

U 

3 

C 

•H 

cr 

■H 

U 

M 

CO 

x: 

■c 

J3 

JS 

o 

CJ 

a 

o 

c 

c 

C 

c 

CO 

CO 

CO 

CO 

>-j 

)-j 

u 

V4 

« 

pq 

PQ 

M 

B 

1-4 

g 

e 

e 

CO 

1 

CO 

1 

0) 

CO 

CO 

1 

0) 

1 

(U 

1 

-a 

"O 

"O 

"D 

•H 

•H 

•H 

•rl 

to 

f 

to 

1 

CO 

to 

1 

0) 

1 

1 

1 

rH 

iH 

rH 

rH 

60 

60 

60 

60 

c 

C 

c 

c 

•W 

•H 

•H 

•H 

Crt 

w 

C/3 

W 

x; 

x: 

,C 

o 

o 

o 

c 

C 

c 

CO 

to 

CO 

Vl 

)-J 

V4 

pq 

pq 

pq 

g 

g 

e 

CO 

to 

CO 

I 

( 

1 

<u 

o 

0) 

■a 

■o 

XI 

rC 

•H 

•H 

•H 

O 

CO 

CO 

CO 

c 

1 

I 

I 

to 

<u 

OJ 

0) 

M 

rH 

iH 

H 

pq 

60 

60 

60 

c 

c 

C 

(U 

•rl 

•H 

•H 

>, 

Crt 

tW 

M 

S 

s 


Csl 


m 


I 


00 


a\ 


42 


^«psr- 


c 

1 

0 

•H 

o 

tJ 

jr. 

OJ 

A 

iJ 

TO 

01 

a 

x: 

•H 

■u 

M 

^ 

U-l 

(0 

o 

14.4 

4J 

• 

flj 

X 

01 

!-i 

Ol. 

M 

Ql 

•H 

3 

> 

U 

W 

« 

a 

O 

C 

o 

m 

Q 

w 

•H 

>> 

iJ 

0) 

x) 

a 

x: 

TO 

(11 

u 

01 

«] 

4J 

)-< 

4J 

cd 

T) 

CO 

a 

g; 

x: 

TO 

H 

'J 

c 

01 

cd 

X! 

• 

M 

3 

M 

^ 

4J 

C 

0 

0 

6.C 

•H 

O 

C 

4J 

CO 

■H 

V 

W 

0! 

0) 

CJ 

CO 

-C 

01 

~-t 

C 

u 

c 

M 

o 

a 

• 

a 

4J 

^ 

u 

en 

5 

o 

V) 

o 

—i 

TO 

-w 

^ 

rH 

-n 

tn 

SO 

OJ 

a 

01 

X. 

r-i 

r^ 

»J 

60 
•H 

w 

»M 

X 

o 

o 

1) 

u 

x: 

P^ 

TJ 

&. 

01 

« 

T3 

c 

M 

0) 

•H 

60  .H 

O 

O 

r-H 

■>-) 

■u 

•H 

o 

U 

rj 

jz 

■a 

3 

(u 

O 

+-I  x: 

o 

CO 

iTi 

g 

CO 

H 

(U 

o 

>J  l*^  ^ 

3 

O. 

60 

01 

TO 

•H  x; 

M 

fe 

u 

60 

43 


approximately  1  cm  on  center  with  the  upstream  branch  stem  to 
accommodate  the  glass  connecting  tube  with  a  snug  and  flush  internal 
fit.   Duco  cement  was  used  as  a  leak  seal  around  the  connecting  tube 
on  the  upstream  block  faces.   Figure  5  shows  the  connecting  tube  joined 
to  test  section  6.   Branch  outlets  were  mounted  in  a  similar  fashion,  and 
very  satisfactory  results  were  obtained  using  Plexiglas  tubing  for  these. 


C.   The  Optical  Syst 


em 


There  are  three  basic  requirements  to  be  considered  in  designing 
the  optical  system.   First,  this  system  must  be  oriented  such  that  the 
light  path  is  perpendicular  to  the  direction  of  flow  being  studied. 
Secondly,  it  should  be  possible  to  produce  incident  circularly  polarized 
light  as  well  as  plane-polarized  light  of  arbitrary  angular  orientation. 
Finally,  it  should  be  possible  to  photograph  the  resulting  birefrigence 
patterns.   Figure  3  illustrates  the  optical  scheme  adopted. 

The  light-diffusing  screen,  consisting  of  a  60  cm  square  plate 
of  smoked  glass,  was  mounted  horizontally  on  runners  between  the  two 
laboratory  tables.   Mounted  in  this  fashion,  the  screen  provided  both  a 
uniformly  lighted  background  and  a  rigid  level  support  for  the  polarizer 
and  test  sections. 

Two  Sylvania  flood  lamps  (type  R-32,  color  temperature  3200°K, 
375  watt),  constituting  a  white^ight  source,  were  mounted  on  ring  stands 
which  rested  on  the  floor.   These  lamps  could  be  easily  raised  to  compen- 
sate for  intensity  loss  due  to  aging,  and  heat  could  be  dissipated  by  a 
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fan.   A  white-light  source  was  chosen,  because  it  was  felt  that  the 
resulting  colored  bands  would  be  easier  to  distinguish  in  the  final 
photographs  than  the  light  and  dark  bands  obtained  with  monochromatic 
light. 

As  Figure  3  shows,  the  polarizer  was  mounted  directly  on  the 
light-diffusing  screen.  Plane-polarized  light  of  various  orientations 
was  produced  using  the  edge  of  the  screen  as  an  azimuthal  angular 
reference.   A  second  polarizing  filter  mounted  directly  on  the  camera 
lens  served  as  the  analyzer. 

Circularly  polarized  light  could  be  produced  by  inserting 
quarter-wave  plates  (not  shown)  between  the  polarizer  and  test  section 
and  between  the  test  section  and  analyzer  respectively.   Again,  the  edge 
of  the  diffusing  screen  was  used  as  an  angular  reference  for  properly 
inserting  these. 

The  camera  used  was  a  Bolex  Paillard  H'^ie  (16  mm)  Reflex  motion- 
picture  camera  equipped  with  a  25  mm  lens  (1  :  1.4  Switar  Rx) .   This  lens 
in  conjunction  with  5  mm  and  10  mm  extension  tubes  permitted  photography 
from  a  distance  of  a  few  centimeters.  The  film  used  was  Ektachrome  (ER 
type  B)  and  a  Weston  light  meter  (Model  745)  was  used  to  check  light 
intensity  at  the  camera.   The  camera  was  mounted  on  a  tripod  taped  securely 
to  the  floor.   This  arrangement  permitted  easy  focusing  using  the  tripod 
crank.   Once  focused,  the  camera  could  be  locked  in  position  and  the 
reflex  viewing  used  to  accurately  check  on  proper  analyzer  orientation 
by  light  extinction  under  zero-flow  conditions. 


SECTION  IV 

EXPERIMENTAL  TECHNIQUE 

The  experimental  work  was  composed  of  four  parts:   A. 
Preparation  and  storage  of  the  dye  suspension,   B.   Determination  of 
fluid  physical  properties,   C.   Flow  regulation  and  collection  and  D. 
Photography. 

A.   Preparation  and  Storapie  of  the  Dye  Suspension 

The  doubly  refracting  liquid  used  was  a  1.4%  aqueous  suspension 
of  the  aniline  dye  Milling  Yellow  (obtained  from  National  Aniline 
Division,  Allied  Chemical  Company).   The  suspension  was  prepared  directly 
from  the  commercial  product,  a  yellow  powder,  using  the  method  suggested 
by  Peebles,  Garber  and  Jury.^^ 

Distilled  water,  about  50%  in  excess  of  that  desired  for  the  final 
volume,  is  heating  to  boiling.   To  this  is  added  a  thick  paste  of  dye 
powder  and  distilled  water,  sufficient  to  make  a  1%  suspension  by  weight. 
Boiling  and  agitation  are  allowed  to  proceed  until  excess  water  has  evapo- 
rated.  Since  birefringent  sensitivity  is  markedly  affected  by  dye  concen- 
tration, the  final  suspension  should  contain  1.3  to  1.6%  dye  for  best 
results.   The  double^^refraction  property  may  be  tested  during  evaporation 
by  withdrawing  a  small  sample  into  a  test  tube,  cooling  to  room  temperature 
and  observing  whether  or  not  interference  patterns  are  produced  when  the 
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tube  is  agitated  between  illuminated  crossed  polarizing  plates. 

After  the  dye  suspension  was  prepared  in  sufficient  quantity, 
it  was  stored  in  sealed  polyethylene  tanks  until  needed  for  use.   It 
was  found  that  suspensions  stored  in  this  fashion  were  very  stable  with 
respect  to  particle  agglutination  over  a  period  of  two  weeks,  after 
which  some  coagulation  occurred. 

B.   Determination  of  Fluid  Physical  Properties 

To  date,  the  most  extensive  physical-property  studies  of  Milling- 
Yellow  suspensions  are  those  of  Peebles,  Prados  and  Honeycutt^^  estab- 
lishing the  birefringent  and  rheological  behavior.   Their  studies  have 
shown  that  viscosity  data  for  Milling-Yellow  suspensions  are  accurately 
correlated  by  the  Powell-Eyring  equation  for  pseudo-plastic  materials. 
Little  work  has  been  done  to  establish  other  physical  properties  for 
such  suspensions. 

The  physical  properties  required  for  the  present  study  were 
viscosity,  density  and  static  refractive  index  of  a  1.4%  dye  suspension 
at  25''C. 

Viscosity  data  were  obtained  by  using  the  96  cm  length  of 
precision-bore  3  mm  glass  tubing  as  a  capillary  viscometer.   The  dye 
suspension  was  permitted  to  flow  through  the  tube  under  steady-flow  con- 
ditions for  several  different  but  known  pressure  gradients.   The  corres- 
ponding flow  rates  were  then  measured  and  the  Hagen-Poiseuille  equation 
used  to  calculate  apparent  viscosity  (that  is,  the  viscosity  which  a  New- 
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tonian  fluid  would  have  flowing  through  the  tube  at  the  same  flow 
rate  under  the  same  pressure  gradient) .   Typical  experimental  results 
are  listed  in  Table  3,  Appendix  B,  and  Figure  6  shows  apparent  viscosi- 
ties and  the  corresponding  Reynolds  numbers  plotted  as  functions  of 
flow  rate  Q  in  the  3  mm  diameter  tube.   These  data  are  in  strong  agree- 
ment  with  the  data  reported  by  Peebles,  Prados  and  Honeycutt. 

Fluid  density  was  determined  by  weighing  25  ml  aliquots  of  the 
dye  suspension  at  25°C.   The  density  was  found  to  be  1.00  gm/cm  . 

Static— refractive^index  measurements  were  made  using  a  Bausch 
and  Lomb  ref ractometer.   The  static  refractive  index  at  25''C  was  found 
to  be  1.338  for  the  1.4%  suspension.   Static  refractive  index ,  however , 
should  not  be  confused  with  the  indices  of  refraction  associated  with 
the  ordinary  and  extraordinary  rays  when  the  fluid  is  in  motion.   It  is 
not  clear  just  what  relationship  the  static  index  bears  to  these  two, 
though  Peterlin  and  Stuarf^°  have  taken  it  to  be  the  mean  value  of  the 
dynamic  indices. 


C.   Flow  Regulation  and  Collection 


The  procedure  for  steady-flow  experiments  was  as  follows:   Prior 
to  filling  the  pressure  tank  with  fluid,  short  rubber  tubes  were  fitted 
over  each  flow  outlet  and  clamped  to  prevent  flow.   During  the  experi- 
ments, these  tubes  were  undamped  and  served  to  divert  outlet  flows  to 
large  collecting  beakers  without  altering  flow  distributions  in  the 
test  sections. 
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The  pressure  tank  and  test  section  were  positioned  for 
photography  and  the  tank  filled  to  a  depth  of  20  cm  with  the  Milling- 
Yellow-dye  suspension.   The  laboratory  compressed-air  supply  was  then 
connected  to  the  tank  via  a  three-way  valve  opened  initially  to  the 
atmosphere.   The  pressure  gauge  in  series  between  the  valve  and  tank 
permitted  tank  air  pressure  to  be  monitored  directly. 

For  an  experimental  run,  the  outlet  flows  were  opened  and  the 
valve  partially  closed  until  the  desired  driving  pressure  was  registered 
on  the  pressure  gauge.   Successive  experimental  runs  were  made  by  incre- 
menting the  driving  pressure  in  steps  of  10  mm  of  mercury.   Before  and 
after  each  run,  the  three-way  valve  was  opened  to  the  atmosphere,  thus 
reducing  excess  tank  pressure  and  lowering  the  outlet  flows.   This  also 
permitted  the  liquid  level  in  the  tank  to  be  returned  to  the  20  cm 
level.   Depending  upon  the  flow  rate,  experimental  runs  lasted  for  30 
seconds  or  1  minute.   However,  the  tank  was  of  sufficient  capacitance 
that  the  liquid  level  did  not  drop  more  than  a  few  millimeters  during  any 
experimental  run.   Outlet  flows  were  measured  with  graduated  cylinders. 

For  the  pulsatile-flow  experiments,  the  pressure  tank  was  replaced 
by  the  bellows  pump  described  previously.   Since  the  flow  output  of  this 
pump  was  limited  to  the  flow  wave  already  shown  in  Figure  4,  no  flow 
regulation  was  possible  here. 

D.   Photography 

Prior  to  actually  photographing  the  experimental  birefringence 
patterns  it  was  necessary  to  standardize  several  important  photographic 
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variables  such  as  light  intensity,  framing  speeds,  exposures,  depth  of 
field,  the  film  used  and  processing. 

Several  preliminary  runs  were  made  to  test  these  factors  using 
a  color  standard  made  from  strips  of  Kodak  red,  green  and  yellow  Wratten 
filters.   It  was  found  that  Ektachrome  film  (ER  type  B)  with  a  correspon- 
ding light  source  having  a  color  temperature  of  3200''K  gave  very  good 
color  reproduction  when  proper  exposures  were  used.   All  films  were  pro- 
cessed by  the  same  company  so  as  to  minimize  any  differences  due  to 
processing. 

Experimental  birefringence  patterns  were  photographed  with  the 
camera  shooting  down  at  the  test  section  mounted  in  the  horizontal  plane. 
After  the  test  section  was  positioned,  the  desired  exposure  and  depth  of 
field  were  set  and  the  lights  turned  on.   The  camera  was  then  focused 
and  locked  into  position.   The  reflex  viewing  was  used  to  assure  that 
polarizer  and  analyzer  were  crossed  by  checking  for  maximum  light  extinc- 
tion under  zero-flow  conditions.   Light  intensity  was  checked  at  the 
camera  using  a  Weston  light  meter  (Model  745) . 

The  flow  was  started,  and  after  the  desired  flow  conditions  were 
achieved,  the  resulting  birefringence  pattern  was  photographed.   Steady 
flows  were  photographed  at  16  frames/sec  and  pulsatile  flows  at  64 
frames/sec. 

The  pertinent  photographic  data  for  the  flow  photographs  appearing 
in  the  next  section  have  been  listed  in  Table  4,  Appendix  C. 


SECTION  V 


THE  ANALYSIS  OF  FLOW-PATTERN  PHOTOGRAPHS 
AND  EXPERIMENTAL  RESULTS 


In  this  section,  selected  frame  enlargements  of  the  16  mm 
film  are  presented  together  with  experimental  data  obtained  from  flow 
measurements  and  frame-by-frame  analysis  of  the  film. 

The  section  is  divided  into  three  parts.   In  Part  A,  the 
isochromatic-band  pattern  obtained  for  a  fully  developed,  angularly 
symmetric  laminar  flow  through  a  tube  of  circular  cross  section  is 
considered,  and  a  method  is  proposed  for  calibrating  this  pattern 
against  shearing  stress,  including  the  effects  of  curved  boundaries 
and  differences  in  refractive  index  in  the  analysis.   In  Part  B,  steady- 
flow  data  and  isochromatic  photographs  are  presented,  and  birefringence 
theory  together  with  the  calibration  of  Part  A  is  applied  to  estimate 
shear-stress  distributions  for  flows  around  a  cardiac  catheter  and 
through  a  90°  branch.   Experimental  results  are  then  compared  with 
results  from  the  medical  literature  concerning  wall  trauma  and  athero- 
sclerotic disease  at  sites  of  branching  in  arteries.   In  Part  C,  branch- 
entrance-length  data  for  the  pulsatile -flow  experiments  are  presented, 
and  discussion  is  focused  on  the  areas  of  stasis,  boundary-layer  sepa- 
ration and  increased  shear  readily  observed  in  the  films. 
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A.   Isochromatic  Shear-Stress  Calibration 
_for  Laminar  Tube  Flow 

The  problem  of  isochromatic  calibration  in  flow-birefringence 
work  is  that  of  associating  an  appropriate  value  of  shear  stress  with 
any  given  isochromatic  band  in  the  flow  field.   In  general,  this 
requires  that  a  birefringence  experiment  be  run  under  conditions  of 
known  geometry,  shear-stress  distribution  and  isochromatic-band  distri- 
bution.  A  correspondence  or  calibration  may  then  be  drawn  between  the 
shear  field  and  isochromatic  field  and  this  calibration  used  to  determine 
the  shear  field  in  geometries  where  only  the  isochromatic  field  is  known. 

It  will  be  recalled  from  Section  II  that  previous  studies  have 
demonstrated  experimentally  that  fringe  order  N  is  directly  proportional 
to  the  product  of  T(avg),  the  average  shear  stress  along  the  light  path, 
and  £,  the  length  of  this  path  in  the  fluid.   According  to  Equation  (52) 
that  relationship  was  given  as  N  =  M*£T(avg)  where  M*  is  a  characteristic 
of  the  particular  fluid.   Although  M*  is  not  constant  for  all  fluids,  the 
data  of  Prados^l  show  that  this  is  the  case  for  Milling-Yellow  suspensions. 

It  is  evident  from  Equation  (52)  that  for  a  three-dimensional 
flow  system  in  which  the  light-path  length  (fluid  depth)  I  associated 
with  a  given  band  N  is  constant,  the  shear  stress  T(avg)  associated 
with  this  band  will  also  be  constant.   However,  if  this  light-path 
length  varies  along  the  band,  then  the  corresponding  shear  stress  will 
vary  inversely  as  the  path  length,  and  there  will  no  longer  be  only  one 
appropriate  shear  stress  associated  with  the  given  band.   A  reasonable 
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approach  for  this  situation  would  be  to  establish  some  reference  point 
along  the  band  (perhaps  upstream  from  the  position  of  interest) ,  and 
to  refer  the  stress  at  other  positions  on  the  band  to  the  value  at  this 
point.   This,  of  course,  requires  that  the  light-path  length  be  known 
at  the  downstream  position. 

Previous  experimentalists  utilizing  flow  birefringence  were  able 
to  avoid  this  diffuculty  associated  with  test  sections  of  varying  depth 
by  using  essentially  two-dimensional  systems  (such  as  high-aspect-ratio 
channels) . 

As  will  be  seen  in  Part  B,  the  isochromatic  bands  obtained  in 
the  present  study  represent  a  combination  of  the  above  two  cases  ,  with 
the  undisturbed  upstream  state  correspondT.ng  to  constant  light-path 
length  in  the  fluid  and  the  distorted  downstream  state  to  variable 
length.   For  this  reason,  we  shall  establish  the  upstream-f low  condition 
as  a  reference  state  for  branching  flows  and  turn  now  to  the  problem  of 
calibrating  the  isochromatic-band  distribution  obtained  for  laminar  tube 
flow. 

We  consider  a  steady,  incompressible,  fully  developed  laminar 
flow  through  a  cylindrical  tube.   A  typical  isochromatic-band  distribu- 
tion observed  for  such  a  flow  is  that  shown  in  Figure  7  for  test  section 
1.   The  flow  here  is  from  left  to  right  at  a  Reynolds  number  of  50 
(based  on  apparent  viscosity  and  inside  tube  diameter) ,  and  the  incident 
light  is  plane  polarized  perpendicular  to  the  axis  of  the  tube.   It  is 
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seen  that  the  bands  parallel  the  direction  of  flow  and  are  not  displaced 
from  their  original  radial  position  as  they  pass  downstream. 

In  seeking  to  associate  a  characteristic  value  of  shear  stress 
with  each  band,  the  observer  must  remember  two  things:   First,  successive 
bands  seen  in  the  flow  photograph  are  the  net  result  of  viewing  the 
flow  through  a  curved  wall,  and, hence,  through  varying  liquid  depths  in 
the  radial  direction.   Second,  since  curved  boundaries  are  involved,  it 
is  necessary  to  consider  the  possibility  of  optical  distortion  arising 
in  the  flow  photograph  due  to  differences  in  refractive  index  between 
the  different  media  involved. 

A  cross-sectional  representation  of  the  over-all  situation  is 
shown  in  Figure  8.   For  convenience,  the  coordinate  origin  0  has  been 
chosen  at  the  center  of  the  tube  with  lines  0-Ni  and  0-No  representing 
normals  to  the  tube  surface  at  points  B  and  C  respectively.   Each  of 
the  three  media  involved  is  characterized  optically  by  a  static  index 
of  refraction  (n^,  n2,  n3).   The  path  of  a  typical  light  ray  as  it 
enters  the  fluid  from  below  and  passes  to  the  observer  is  indicated  by 
line  DCBA,  and  it  will  be  noted  that  with  the  exception  of  the  ray 
passing  through  the  center  of  the  tube,  all  rays  seen  by  the  observer 
pass  obliquely  from  one  medium  to  the  next,  thereby  undergoing  refraction. 
This  refraction  is  characterized  by  the  angles  G^^,  02,  63  and  0a  which 
are  predictable  from  geometrical  considerations  and  Snell's  law  of 
refraction.   These  relationships  will  be  considered  presently.   The 
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OBSERVER 
A  A  ^A 


MEDIUM  I 
(n,) 


LIGHT- DIFFUSING    SCREEN 


Figure  C.  Schematic  cross  section  of  a  circular  tube,  illustrating 
the  refraction  of  a  typical  light  ray  passing  from  the  fluid  to  the 
observer.  Light  path  DCBA  is  characterized  by  angles  GisG2>  ^3  and 
04.   These  angles  are  determined  by  the  narametsr  K  and'^Snell's  law 


of  refract 
are  norinals 


I,   n2  and  n3  are  indices  of  refraction. 


Nj  and  N2 
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observer  (camera)  viewing  from  above  sees  only  those  light  rays 
emerging  parallel  to  the  y-axis,  and  this  implies  that  chords  such  as 
DC  ,  which  are  the  light  paths  in  the  fluid  corresponding  to  isochro- 
matic  bands  seen  the  the  flow  photograph,  will  not  in  general  be  parallel 
to  the  y-axis. 

For  the  present  analysis,  we  shall  choose  chord  DC  as  being  a 
typical  chord  associated  with  an  isochromatic  band  seen  at  point  A  in 
the  photograph  and  shall  identify  this  chord  by  its  perpendicular  dis- 
placement from  the  origin.   The  lateral  displacement  from  the  y_axis 
of  the  band  at  point  A  is,  of  course,  dependent  on  the  chord's  displace- 
ment from  the  origin,  and  this  dependence  will  now  be  established  by 
considering  the  above-mentioned  relationships. 

It  should  be  mentioned  that  there  is  some  question  as  to  the 
justification  for  using  Snell's  law  of  refraction  for  a  doubly  refrac- 
ting material,  since  it  is  kno\<m  that  the  ordinary  ray  obeys  this  law 
but  the  extraordinary  ray  does  not.   In  addition,  each  of  these  rays 
has  its  own  dynamic  index  of  refraction,  so  the  use  of  the  static  index 
must  also  be  justified.   Actually,  these  questions  pose  no  difficulty 
here,  because  both  rays  travel  the  same  path,  and,  as  Peterlin  and 
Stuart  have  pointed  out,    the  difference  between  the  three  indices  is 
very  small.   Therefore,  for  the  purpose  of  determining  the  light  path, 
the  static  index  and  Snell's  law  may  be  used. 

Referring  again  to  Figure  8,  one  may  write  the  following  relations 
between  angles  0-j^,  02,  O3  and  0^; 


58 


ni  sin  0^  =  n2  sin  02       by  Snell's  law  (53) 

n2  sin  03  =  n3  sin  04       by  Snell's  law  (54) 

R  cos (03-02)  +  BC  cos  02  =  KR   Geometry  (55) 

R  sin (02-02)  =  BC  sin  02        Geometry  (55) 

If  0^  is  considered  the  independent  variable  (this  can  be  done  because 
of  the  principle  of  optical  reversibility^^),  then  these  equations  may 
be  solved  simultaneously  to  give 

-1  "1 

02  =  sin  (n^  sin  0^)  (57) 

-1   °1  ,• 

03  =  sin   (K^psi"  01)  (58) 

-1   "1 

04  =  sin  (K—  sin  q^)  (5g) 

-where  the  parameter  K  is  simply  (outside  diameter) /(inside  diameter). 
For  the  special  case  of  a  hole  drilled  parallel  to  a  flat  face  of  a 
Tlexiglas  block,  these  equations  still  apply  providing  K  is  given  a 
numerical  value  of  unity  and  n-,^  is  replaced  by  n  .   This  is  equivalent 
to  immersing  the  observer  in  the  Elexiglas,  since  perpendicular  viewing 
would  create  no  refraction  at  the  air-block  interface. 

One  notes  from  Equation  (59)  that  the  angle   04  is  completely 
independent  of  medium  2  (that  is  n2) .   This  is  a  well-known  result  in 
optics.   Another  interesting  and  important  observation  is  that  the 
diameter  ratio  K  plays  a  central  role  in  determining  how  much  of  the 
flow  field  is  actually  visible  to  the  observer.   This  fact  may  be 
appreciated  by  considering  the  following  relationships,  which  may  be 
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deduced  from  equations  (57)  through  (59)  and  the  fact  that  in  this 
study  n2>no>n-, . 


IF 


K  =  1. 0  and  n^  =  n2 
(Plexiglas-block  case) 


n- 


1  <K  <- 


ni 


K  = 


n. 


THEN 


e^  >  01  =  03  =  02 


0^  >  0^   >  03  >  02 


0^  =  0,  >  03  >  ©2 


(60) 


(61) 


(62) 


"3 


<  K  < 


^1        ^1 


04  >  01  >  03  >  02       (63) 


K  = 


04  >  01  =  03  >  02 


(64) 


K  > 


04  >  ©3  >  01  >  02       C65) 


These  relationships  show  that  with  the  exception  of  the  second  case  (61), 
0^  is  always  greater  than  or  equal  to  0i.   This  is  important,  because  as 
shown  in  Figure  8,  the  magnitude  of  0^  must  range  from  0  to  90°  if  the 
observer  is  to  see  the  entire  flow  field.   At  the  same  time,  it  will  be 
noted  that  ©i  physically  can  be  no  greater  than  90°.   If  case  (61)  holds, 
that  is,  if  K  falls  between  1.0  and  n3/ni,  then  04  will  be  less  than  Q-.^ 
and  the  observer  will  not  be  able  to  see  the  entire  flow  field  but  only 


-1 


-1.  ^ 


the  part  of  it  for  which  0,  ranges  from  -  sin   (K — —)    to  +  sin""  (K ~) 

as  dictated  by  Equation  (59).   Thus,  an  experimenter  conducting 
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flow-birefringence  experiments  employing  circular  tubes  and  wishing 
to  visualize  the  entire  flow  field  must  consider  the  geometrical  range 
1  <  K  <  n^/n^^  a  forbidden  zone  when  choosing  experimental  test 


sections, 


For  the  flow  pictured  in  Figure  7,  the  appropriate  numerical 
values  are  K  =  3.0,  n^  =  1.00,  n^   =  1.5  and  n^  =   1.338.   Since   ■ 
this  value  of  K  is  greater  than  n3/n3^,  one  is  assured  that  the  entire 
flow  field  is  visible  in  that  photograph.   In  fact,  the  entire  flow 
field  is  visible  in  all  photographs  appearing  in  this  study.   For  the 
branching-flow  figures  the  appropriate  numerical  values  are  K  =1.0, 
"l  "  "2  "  ^-55  ^^^  ^3   =  1.338. 

It  can  be  seen  from  Figure  8  that  the  perpendicular  displacement 
of  the  typical  chord  DC  from  the  origin  is  given  by  R  sin  0^,  and  that 
the  lateral  displacement  of  the  corresponding  band  at  point  A  is  KR 
sin  Q-^.      Defining  a  distortion  ratio   as  (band  displacement) /(chord  dis- 
placement), one  obtains  a  distortion  ratio  of  K  l^^L—i—  .   From 

sin  0^ 

Equation  (59)  the  value  of  sin  0^/  sin0^  is  n3/Kn^.   Substituting  this 
above  gives  a  distortion  ratio  of  n3/n;L-   Thus,  our  analysis  shows  that 
for  circular  geometry  the  distortion  ratio  is  constant  and  independent 
of  both  the  dimensions  and  the  material  of  medium  2  (the  tube). 
Furthermore,  this  distortion  represents  a  magnification  effect  with  the 
distortion  ratio  n^/n^   as  the  magnification  factor.   In  the  present  study, 
this  factor  has  the  value  1.338  for  the  glass  tubes;  but  since  n-^  must  be 
replaced  by  n2  for  the  Plexiglas-block  studies,  the  factor  becomes  0.865 
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for  them.   In  either  case,  the  only  optical  distortion  arising  from 
curvature  is  magnification  or  shrinking  of  the  original  image  by  a  con- 
stant known  factor.   We  can  use  this  information  to  obtain  one  more 
important  result.   Let  the  band  displacement  at  point  A  be  denoted  by 
"b"  and  the  corresponding  chord  displacement  by  "a".   Then,  according 
to  the  above  discussion,  b  =  an-j/n-,.   The  displacement  b  will  have  its 
maximum  value  b(max)  when  a  has  its  maximum  value  a (max)   =  R.   Thus, 
b(max)   =  Rn^/nj^.   If  both  displacements  are  now  made  dimensionless  by 
normalizing  them  to  their  maximum  values  and  the  fraational  displacements 
"3"  and  "a"  defined  as  g  =  b/b(max  and  a  =  a/a (max),  we  obtain 
the  important  result  that  6  =   (an^/n^^)/ (Rn3/n, )   =  a/R  =  a.   This 
means  that  the  fractional  displacement  (dimensionless  radius)  of  an 
isochromatic  band  measured  in  the  flow  photograph  is  precisely  equal  to 
that  of  the  corresponding  chord  (light  path)  measured  inside  the  tube. 

Now,  having  related  the  geometrical  properties  of  the  band  in  the 
flow  photograph  to  those  of  the  corresponding  light  path  in  the  fluid, 
we  are  in  a  position  to  derive  the  expression  relating  the  fractional 
displacement  a  (or  3)  to  the  average  shear  stress  along  the  chord  DC. 

Up  to  this  point,  no  mention  has  been  made  of  the  nature  of  the 
shear„stress  distribution  in  the  fluid.   From  here  on,  however,  it  will 
be  assumed  that  the  flow  field  in  the  tube  is  angularly  symmetric,  which 
implies  that  the  stress  field  is  also  angularly  symmetric.   The  assumption 
of  angular  symmetry  was  verified  experimentally  for  the  band  distribution 
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shown  in  Figure  7. 

In  Section  II  it  was  shown  that  for  the  steady,  incompressible, 
angularly  symmetric  and  fully  developed  flow  being  considered  here 
the  equations  of  motion  could  be  integrated  to  give  the  shear-stress 
field  as  a  linear  function  of  radial  position  r.   That  function  was 
given  by 

X  =  v^  ,  (66) 

8P    R 
where  the  wall  stress  x^  is  -  (-5 — )  —j—  . 

The  problem  of  determining  the  average  value  of  x  over  the  light 
path  DC  can  be  considerably  simplified  by  taking  advantage  of  the 
— angular  symmetry. 

- -Consider  the  enlarged  view  of  the  flow  cross  section  shown 
schematically  in  Figure  9.   Though  the  light  path  DC  might  in  general 
occupy  a  position  such  as  that  indicated  by  the  position  vector  r,  the 
assumption  of  angular  symmetry  permits  this  chord  to  be  rotated  about 
the  origin  to  an  arbitrary  new  position  without  in  any  way  altering 
the  shear-stress  distribution  along  the  chord.   In  particular,  if  one 
chooses  to  rotate  chord  DC  to  the  position  D'C'  parallel  to  the  y_axis, 
then  the  average  stress  from  D'  to  C  is  the  same  as  the  average  stress 
from  y  =•  0  to  y  =  y*,  owing  to  symmetry  with  respect  to  reflection 

in  the  x-^axis.   Thus,  the  average  shear  stress  along  the  chord  can  be 

2  2  % 
obtained  by  integrating   x  from  y  =  0  to  y  =  y*  =  (R  -a  )  ,  where 


:^?-W^''  ■■' 


^.U-'^'^' 
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X 


Figure  9.   Schematic  flow  cross  section,  illustrating  a  typical  light 
path  along  which  the  average  shear  stress  is  required.   For  angularly 
symmetric  flow,  the  average  stress  along  light  path  DC  is  ecual  to 
that  along  path  D'C'. 


"■?^ 
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the  parameter  a  is  the  perpendicular  displacement  of  the  chord  from 
the  origin. 

Along  the  path  D'C',  Equation  (66)  takes  the  form 

^  =  ^w 1 (67) 

and  the  average  shear  stress   T(avg)  is  given  by 


Tdy 
T(avg)  =  rL^ ,  y*  =   (R-a)^  (68) 


dy 


Substituting  x  from  Equation  (67)  into  Equation  (68)  and 
performing  the  indicated  integrations  gives 


"w   (        a^        ^   ^  (R^-a^^ 

t(avg)   =  -^^  ^  1  + ^       In  [ 3  }   (69) 

R(r'^  -  a^)' 


+  R     ) 


If  we  now  define  the  dimensionless  variable  $  as  T(avg)  /x 

and  rewrite  Equation  (69)  in  terms  of  the  fractional  displacement  a, 
we  obtain  finally 


■^w        2  J     (l-a^)%        ^     a     J 


(70) 


This  equation  gives  the  desired  average  shear  stress  along  the 
light  path  as  a  fraction  of  the  wall  stress  in  the  tube.   On  the  basis 
of  physical  intuition,  one  would  expect  $  to  converge  to  the  values  of 
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0.5  and  1.0  for  a  =  0  and  a  =  1  respectively.  Although  inspection 
of  Equation  (70)  shows  that  $  is  indeterminate  at  both  endpoints,  a 
single  application  of  L'Hospital's  rule  is  sufficient  to  demonstrate 
that  $  does  indeed  converge  to  the  expected  values. 

The  length  £  of  the  chord  D'C  (or  DC)  may  be  written  in  terms 
of  a  by  considering  the  geometry  of  Figure  9.   The  relationship 
between  £  and  a  is 

I  /2R  =  (1  -a^)^  (71) 

The  right  side  of  Equation  (71)  thus  gives  the  light -path  length  in  the 
fluid  as  a  fraction  of  the  inside  tube  diameter,  and,  as  expected,  this 
length  has  a  value  of  one  diameter  for  a  =  0  and  zero  for  a  =  1. 
Numerical  results  from  Equations  (70)  and  (71)  have  been 
summarized  in  Table  2.   The  first  column  in  the  table  represents  either 
the  fractional  displacement  of  the  light  path  from  the  center  of  the 
tube  or  the  fractional  displacement  of  an  isochromatic  band  as  measured 
in  the  flow  photograph.   The  second  column  gives  the  length  of  the  corres- 
ponding light  path  in  the  fluid  in  diameters,  and  the  third  column  gives 
/  the  average  shear  stress  along  the  light  path  as  a  fraction  of  the  wall 
stress  in  the  tube. 

As  an  example  of  the  use  of  this  table,  consider  the  green  band 
which  was  seen  in  Figure  7.   In  the  photograph,  the  center  of  this  band 
is  displaced  10.4  mm  from  the  center  of  the  tube.   The  edge  of  the 
birefringence  pattern  is  displaced  13.0  mm  from  the  center.   The  fractional 
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Table  2 

LIGHT- PATH  LENGTH  AND  AVERAGE  SHEAR  STRESS 
ALONG  THIS  PATH  AS  DIMENSIONLESS  FUNCTIONS 
OF  FRACTIONAL  DISPLACEMENT  IN  THE  FLUID  MEDIUM* 


Fractional 

Dimensionl 

.ess 

Dlmensionless 

Displacement,  (a  or 

B)   Path  Length, 

(  V2R) 

Average 

Shear  Stress,  $ 

0.000 

1.000 

0.500 

0.050 

0.999 

0.505 

0.100 

0.995 

0.515 

0.150 

0.989 

0.529 

0.200 

0.980 

0.547 

0.250 

0.968 

0.567 

0.300 

0.954 

0.588 

0.350 

0.937 

0.612 

0.400 

0.917 

0.637 

0.450 

0.893 

0.663 

0.500 

0.866 

0.690 

0.550 

0.835 

0.718 

0.600 

0.800 

0.747 

0.650 

0.760 

0.777 

0.700 

0.714 

0.807 

0.750 

0.661 

0.838 

0.800 

0.600 

0.870 

0.850 

0.527 

0.902 

0.900 

0.436 

0.934 

0.950 

0.312 

0.967 

1.000 

0.000 

1.000 

*N0TES:   a.   This  calibration  applies  only  for  fully  developed 

angularly  symmetric  laminar  flow  in  a  circular  tube. 

b.  Fractional  displacement  a  measures  the  perpendicular 
displacement  between  the  center  of  the  tube  and  the 
chord  representing  the  light  path.   Fractional  dis- 
placement &  measured  in  the  flow  photograph  for  the 
corresponding  isochromatic  band  is  exactly  equal  to  c 

c.  £/2R  is  calculated  from  Equation  (71). 

d.  $  is  calculated  from  Equation  (70). 
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displacement  of  this  band  is  therefore  3  =  10.4/13.0  =  0.8.   Consulting 
Table  2.  one  finds  that  the  light-path  length  in  the  fluid  corresponding 
to  this  displacement  is  0.600  diameters.   Since  the  tube  is  known  to 
have  an  inside  diameter  of  3.00  mm.  the  light  path  is  1.80  mm  long. 
Finally,  from  Table  2,  the  average  shear  stress  associated  with  this 
band  is  0.870  times  the  know  wall  stress  of  90  dynes/cm2.   Therefore, 
the  shear  stress  associated  with  the  green  band  is  78.3  dynes/cm2. 

It  should  be  emphasized  that  the  calibration  given  by  Table  2 
applies  only  to  isochromatic  bands  obtained  for  tubes  of  circular- 
cylindrical  geometry  under  the  steady  laminar-flow  conditions  described 
earlier.   These  are  the  upstream  conditions  prevailing  in  the  steady-flow 
experiments  of  this  study. 

-    --The  problem  of  determining  the  shear  stress  at  downstream 
-entrance-flow  positions  in  terms  of  the  upstream  conditions  will  be 
considered  in  Part  B  with  numerical  examples  for  flows  around  a  cardiac 
catheter  and  through  a  90"  branch. 

B.   Steady-Flow  Experiments 

To  facilitate  discussion  of  the  following  figures,  it  would  be 
appropriate  to  discuss  first  the  conventions  which  have  been  followed 
regarding  the  orientation  of  the  flow  photographs  and  general  terminology. 

The  photographs  appearing  in  this  section  show  typical  isochroma- 
tlc^band  patterns  for  flows  through  both  straight  and  branched  test 
sections.   The  observer  should  keep  in  mind  that  for  each  photograph 


68 


the  flow  is  steady  and  from  left  to  right  and  for  branched  test 
sections  this  flow  is  such  as  to  impinge  on  the  bifurcation  (that  is, 
all  branched  flows  are  of  the  diverging  type). 

The  convention  followed  in  labeling  branches  is  that  side  arms 
have  been  designated  "B"  and  downstream  main  stems  "A",  with  the  term 
"upstream"  referring  to  the  pre-branch  segment. 

All  Reynolds  numbers  have  been  based  on  apparent  viscosity  and 
the  pertinent  tube  diameter  as  a  characteristic  length  (for  the  single 
branch  of  rectangular  cross  section,  this  length  has  been  taken  as 
four  times  the  hydraulic  radius,  defined  as  the  cross-sectional  area 
of  the  stream  divided  by  the  wetted  perimeter) . 

The  term  "flow  ratio"  applies  only  to  branched  test  sections  and 
is  defined  as  the  ratio  of  the  volumetric  flow  rate  in  the  side  arm  to 
that  for  the  downstream  main  stem  (i.e.,  Qg/Q^) •   Similarly,  the  term 
"area  ratio"  is  defined  as  the  ratio  of  the  cross-sectional  area  of 
the  stream  in  the  side  arm  to  that  for  the  downstream  main  stem  Ci.e., 
Ag/A^) . 

Unless  otherwise  noted,  the  values  of  wall  stress  and  Reynolds 
number  given  in  the  figures  refer  to  the  upstream  (pre-branch)  values 
of  these  quantities.   In  each  case,  upstream  wall  stress  has  been 
experimentally  determined  as  the  negative  of  the  pressure  gradient  mul- 
tiplied  by  half  the  inside  radius  at  the  upstream  position. 

Finally,  all  entrance-length  measurements  have  been  expressed  in 


J' 
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dimensionless  form  by  normalizing  them  to  the  diameter  of  the  branch 
for  which  they  apply  (i.e.,  side-arm  entrance  lengths  have  been 
normalized  with  respect  to  side-arm  diameter  and  main-stem  entrance 
lengths  with  respect  to  main-stem  diameter). 

Prior  to  the  actual  branching  studies,  it  was  decided  to  employ 
the  birefringence  technique  to  obtain  isochromatic  patterns  for  steady 
flows  around  various  needles  and  cardiac  catheters  inserted  radially 
and  axially  into  a  circular  tube.   Since  the  birefringence  technique 
had  not  been  used  previously  for  this  purpose,  it  was  hoped  that  useful 
information  would  be  obtained  regarding  the  creation  and  propagation 
of  flow  disturbances  caused  by  these  devices. 

Figure  10  shows  an  isochromatic-band  distribution  for  steady 
flow  around  a  catheter  (O.D.  =  1.67  mm)  inserted  axially  into  a  cylin- 
drical tube  (I.D.  =  3.00  mm).   Here,  plare-polarized  light  has  been 
used  with  the  plane  of  polarization  perpendicular  to  the  axis  of  the 
tube.   The  upstream  Reynolds  number  for  the  experiment  was  10,  and  the 
wall  stress  was  75  dynes/cm^. 

In  the  upstream  region,  the  bands  are  distributed  much  as  they 
were  for  the  fully  developed  flow  in  Figure  7,  but  as  they  encounter  the 
catheter,  these  bands  become  noticeably  displaced  toward  the  tube  wall. 
The  flow-interference  front  caused  by  the  catheter  extends  approximately 
one  diameter  upstream,  and  the  region  between  the  catheter  and  the  tube 
wall  is  recognized  as  one  of  accelerated  flow  owing  to  the  restriction  of 
flow  area.   The  band  crowding  seen  in  this  region  is  typical  of  the 
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patterns  seen  for  regions  of  accelerated  flow  and  will  also  be  seen 
in  photographs  of  branching  flows.   This  crowding  of  bands  provides  a 
clue  to  areas  where  one  would  expect  to  find  increased  velocity 
gradients  and  shearing  stress. 

It  is  apparent  from  Figure  10  that  the  flow  in  the  entrance 
region  surrounding  the  catheter  does  not  correspond  to  the  fully 
developed  flow  discussed  in  Part  A.   This  means  that  the  shear-stress 
calibration  given  by  Equation  (70)  in  Part  A  cannot  be  applied  directly 
to  calculate  the  stress  distribution  for  the  entrance  region.   In  fact, 
since  the  light-path  length  varies  along  the  bands  in  this  region,  there 
is  no  unique  stress  which  can  be  associated  with  any  band.   However,  it 
was  intimated  in  Part  A  that  the  fundamental  Equation  N  =  M*JlT(avg) 
could  be  used  to  relate  the  stress  distribution  in  the  entrance  region 
to  that  in  the  upstream  region  of  fully  developed  flow,  where  the  cali- 
bration of  Part  A  does  apply.   In  particular,  since  the  coefficient  M* 
has  been  found  to  be  independent  of  shear  rate  for  Milling-Yellow  sus- 
pensions, one  can  apply  the  above  equation  to  an  upstream  and  downstream 
point  on  the  same  band  (or  on  two  different  bands  of  the  same  order  N) 
to  obtain 

T(avg.  down)  =   [  ^^ll^^      T(avg,  up)  (72) 

From  Equation   (72)  it  is  evident  that  the  average  shear  stress 
along  the  light  path  at  the  downstream  position  can  be  calculated  from 
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the  upstream  value  for  the  given  band  providing  that  the  light-path 
length  at  each  of  the  txjo  points  is  known. 

It  was  shown  in  Part  A  that  for  flow  through  tubes  the  light-path 
length  Z   is  related  directly  to  the  observed  dimensionless  displacement 

e  of  the  isochromatic  band  at  the  point  in  question  by  Equation  (71), 

2  % 
Z/2R   =  (1-e  )  ,  where  we  have  used  the  fact  that  B  =  a.   Since  this 

equation  was  derived  solely  from  geometrical  considerations,  it  may  be 

applied  accurately  to  both  upstream  and  downstream  positions  and 

Equation  (72)  written  as 


(l-B^(up))' 
(l-B^(down))' 


T(avg,  dora)  =— —^   T(avg,  up)  (73) 


— Here  6(up)  and  6(down)  are  the  dimensionless  displacements  of  the 
isochromatic  band  measured  in  the  flow  photograph  at  the  upstream  and 
downstream  points  respectively. 

Rewriting  Equation  (73)  in  terms  of  the  upstream  wall  stress 
\,(up)  and  the  dimensionless  average  shear  stress  ^(up)  one  has 


2.   ..% 

(74) 


Z.  2 

T(avg,  down)  =   (1-g  (up))     ^^^^^  ^ 
(1-3  (down))'        "" 

All  the  quantities  on  the  right  side  of  Equation  (74)  are  known, 
either  from  the  birefringence  experiment  or  from  the  shear-stress  cali- 
bration of  Part  A.  The  quantities  $(up),  (l-B^(up))'  and  (l-B^(down))^ 
can  be  obtained  either  from  Table  2  in  Part  A  for  the  appropriate  value 
of  B  or  by  direct  calculation. 


-  "^s^: 
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As  an  example,  consider  the  dark  band  appearing  upstream  in 
the  center  of  the  tube  shown  in  Figure  10.   Encountering  the  catheter, 
this  band  is  displaced  from  its  original  position,  B  =  0,  through 
intermediate  values  of  3  to  a  final  position,  6  =  0.5,  between  the  tip 
of  the  catheter  and  the  tube  wall.   Choosing  this  point  as  the  downstream 
point  of  interest  and  consulting  Table  2,  one  finds  the  following: 

6(up)  =  0,  (l-B^(up))^  =  1.000  and  (|>(up)  =  0.500 

B(down)  =  0.5,  (1-3^ (down) ) ^  =  0.866 

The  known  wall  stress  for  this  experiment  is  given  in  Figure  10  as 
75  dynes/cm  .   Inserting  these  values  in  Equation  (74)  gives 

T(avg,  down)  =  (1.000/0.866)  (0.500)  (75)  =  43.3  dynes/cm^. 

For  comparison,  the  average  shear  stress  corresponding  to  this  band  in 
the  upstream  position  is 


T(avg,  up)  =  $(up)  T^(up)   =   (0.500)  (75)   =  37.5  dynes/ 


2 

cm  , 


As  expected,  the  region  of  accelerated  flow  is  also  one  of  higher  shear. 

It  should  be  emphasized  that  the  above  calculation  gives  only 
the  average  shear  stress  along  the  light  path  at  a  single  dc^rastream  . 
position  and  applies  only  when  the  reference  (upstream)  flow  is  fully 
developed  as  described  earlier.   By  repeating  this  calculation  for  bands 
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falling  close  to  the  wall  of  the  tube  one  can  obtain  an  estimate  of 
the  shear  stress  acting  on  the  wall. 

Figure  11  shows  another  experiment  in  which  a  beveled  needle 
(shank  diameter  =  0.92  mm)  has  been  inserted  through  the  wall  of  the 
cylindrical  tube  (I.D.  =  3.00  mm)  and  oriented  with  the  bevel  facing 

upstream.   The  Reynolds  number  for  this  experiment  was  40,  and  the 

2 

upstream  wall  stress  was  80  dynes/cm  .   Again,  plane-polarized  light 

has  been  used  with  the  plane  of  polarization  perpendicular  to  the 
axis  of  the  tube. 

The  photograph  has  been  underexposed  to  emphasize  the  two 
distinct  vortex  patterns  seen  trailing  behind  the  needle.   These 
patterns  persist  for  approximately  two  diameters  downstream  where  the 
disturbance  becomes  damped  out  because  of  fluid  viscosity.   The  region 
of  disturbed  flow  extends  less  than  one  diameter  upstream. 

Though  the  vortex  patterns  are  similar  in  size  and  shape  in 
this  photograph,  motion  pictures  reveal  that  when  the  needle  is  rota- 
ted 90°,  the  lower  vortex  disappears  almost  entirely  while  the  upper 
one  remains  essentially  unchanged.   This  is  evidence  that  the  edge  of 
the  bevel  creates  far  less  disturbance  than  does  the  face.   T-Then  the 
needle  is  rotated  further  until  the  bevel  faces  downstream,  the  lower 
vortex  reappears  but  is  not  as  prominent  as  the  upper  one.   The  same 
experiment  using  a  cylindrical  needle  having  no  bevel  produces  only  a 
single  vortex  pattern  which  appears  in  the  center  of  the  tube  and  is 
much  more  streamlined  than  those  in  Figure  11.   Ttiis  disturbance  also 
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persists  downstream  about  two  diameters. 

The  most  significant  finding  from  these  experiments  was  that 
when  the  needles  were  withdrawn  to  a  position  flush  with  the  wall  of 
the  tube  no  disturbances  could  be  seen  in  the  birefringence  patterns. 
This  strongly  suggests  that  lateral  pressure  measurements  could  be 
made  in  small-diameter  models  as  well  as.  in  vivo   without  appreciably 
disturbing  the  flow.  _ 

Figure  12  shows  an  isochromatic-band  distribution  obtained 
using  circularly  polarized  light  for  the  first  of  the  branching-flow 
experiments.   The  branch  shown  (test  section  6)  has  a  branching  angle 
of  30°  and  an  area  ratio  of  unity  (the  side  arm  and  main  do\>mstream 
stem  both  having  an  inside  diameter  of  3.17  mm).   The  upstream 
Reynolds  number  for  this  experiment  was  116  with  a  corresponding  v/all 
stress  of  110  dynes/cm^. 

Although  the  bands  are  parallel  to  the  axis  of  the  tube  in  the 
upstream  segment,  they  tend  to  converge  toward  the  inner  walls  and 
distribute  themselves  concavely  about  the  bifurcation  point  in  the 
entrance-flow  region.   This  region  of  disturbed  flow  extends  doxmstream 
for  a  few  diameters  in  both  arms  until  the  bands  return  to  the  parallel 
disturbution  indicative  of  laminar  flow.   The  fact  that  the  bands  tend 
to  be  distributed  concavely  about  a  stagnation  point,  such  as  the 
point  of  bifurcation,  is  helpful  in  locating  zones  of  stasis  at  other 
positions  in  a  flow  test  section.   This  will  be  more  apparent  in  later 
photographs.   The  displacement  of  bands  toward  the  inner  walls  of  the 
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branch  indicates  that  the  regions  directly  adjacent  to  the  bifurcation 
represent  areas  of  accelerated  flow  and  increased  shearing  stress, 
while  the  absence  of  bands  toward  the  outer  walls  indicates  boundary- 
layer  separation  in  those  areas.   Aside  from  the  regions  immediately 
adjacent  to  the  bifurcation,  there  is  no  evidence  for  increased  stress 
at  any  other  site  in  this  branch. 

Figure  13  gives  a  comparison  of  the  flow-distribution 
characteristics  for  three  different  30"  branches.   Here,  flow  ratio 
(Qg/Q^)  has  been  plotted  against  upstream  Reynolds  number  with  area 
ratio  (Ag/A.)  as  a  parameter.   The  upper  curve  applies  to  the  branch 
just  discussed,  and  the  lower  two  apply  to  branches  having  side-arm 
diameters  of  2  mm  and  1  mm  respectively. 

From  the  over-aT.1  decrease  of  flow  ratio  with  increasing 
Reynolds  number  shown  by  all  three  branches,  it  is  evident  that  the 
upstream  flow  increasingly  favors  a  straight  path  into  the  main  down- 
stream stem  A  as  the  upstream  fluid  momentum  per  unit  voltmie  increases. 
The  marked  separation  between  successive  pairs  of  curves  in  Figure  13 
simply  verifies  that  reducing  the  side-arm  diameter  is  an  effective 
means  of  controlling  side-arm  flow  over  the  entire  range  of  upstream 
Reynolds  numbers.   It  will  be  seen  later  that  this  control  is  much  better 
than  that  achieved  by  increasing  the  branching  angle.   The  disproportionate 
increase  of  flow  ratio  which  occurs  below  the  Reynolds  number  of  25  is 
londoubtedly  a  manifestation  of  wall  attachment  by  the  flow,  similar  to 
that  observed  by  Fox  and  Hugh-'-^  in  their  experiments  with  open-channel 
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flows.   The  data  of  Figure  13  indicate  that  for  30°  branches  wall 
attachment  can  cause  the  side-arm  flow  to  greatly  exceed  the  main- 
stem   flow,  but  only  if  the  branch  area  ratio  approaches  unity  and 
only  below  a  certain  upstream  Reynolds  number. 

The  phenomenon  of  wall  attachment  (sometimes  called  the  Coanda 
effect)  is  known  to  arise  because  of  a  low-pressure  region  generated 
along  the  outer  side-arm  wall  owing  to  fluid  moving  along  curved 
streamlines  concave  toward  this  wall.   In  branches,  this  low-pressure 
region  tends  to  divert  a  greater  fraction  of  the  incoming  flow  into 
the  side  arm  until  the  upstream  momentum  forces  can  overcome  the  pressure 
forces  causing  the  attachment.   At  still  higher  flow  rates,  the  boundary 
layer  separates  from  the  outer  wall  and  a  vortex  is  generated  there. 
The  data  of  Knox   indicate  that  there  is  a  sharp  drop  in  pressure 
associated  with  this  vortex  and  a  corresponding  rise  of  pressure  in  the 
main  stem  due  to  the  exchange  of  incoming  fluid  kinetic  energy  for  pres- 
sure energy  there. 

It  has  been  found  that  very  reliable  fluid  amplifiers  and  fluidic 
digitalr-logic  devices  can  be  produced  by  carefully  designing  branches 
with  control  ports  at  the  vortex  regions.   Fluid  can  be  introduced 
through  these  ports  to  destroy  a  given  low-pressure  zone,  thus  causing 
the  main  flov;  to  switch  from  one  arm  of  the  branch  to  the  other.   Unfor- 
tunately, very  little  has  been  published  concerning  the  branching-flow 
characteristics  of  fluidics  devices,  since  most  of  this  information  has 
been  considered  proprietary. 
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It  was  decided  to  investigate  further  the  flow-distribution 
characteristics  of  the  30°  branch  shown  in  Figure  12  by  repeating  the 
above  experiment  with  pure  water  as  the  working  fluid  and  to  compare 
these  results  with  those  for  the  equivalent  experiment  with  Milling- 
Yellow  dye  (the  upper  curve  of  Figure  13) . 

Figure  14  shows  the  experimental  results  obtained.   Three  things 
are  immediately  apparent  from  this  figure.   First,  both  the  dye  and 
water  show  an  inverse  dependence  of  flow  ratio  on  upstream  Reynolds 
number.   Second,  though  these  curves  have  the  same  general  shape  and 
apply  to  the  same  branch,  they  are  separated  by  a  decided  jump.   Third, 
though  wall  attachment  causes  the  dye  stream  to  exceed  a  flow  ratio  of 
unity  only  for  Reynolds  numbers  below  25,  the  same  effect  occurs  in 
water  at  much  higher  Reynolds  numbers;  and, in  fact,  below  a  Reynolds 
number  of  250  in  water  the  wall -attachment  forces  are  sufficiently  strong 
so  that  no  flow  enters  the  main  downstream  stem  (that  is,  complete  main- 
stem  cutoff  occurs) . 

These  differences  can  be  best  understood  by  comparing  the  two 
fluids  on  the  basis  of  upstream  momentum  per  unit  volume  rather  than  on 
the  basis  of  upstream  Reynolds  number.   Since  the  Reynolds  number  reflects 
viscous  forces  as  well  as  momentum  forces,  the  reduction  in  dye  viscosity 
with  increasing  flow  rate  affects  the  lower  flow-ratio  curve  in  Figure 
14  while  no  such  effect  occurs  with  water.   Recalling  that  over  the  flow 
range  covered  by  these  experiments  the  dye  suspension  had  a  viscosity 
from  ten  to  fifteen  times  that  of  the  water,  then  when  the  Rejmolds 
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Figure  14.   Flow  ratio  versus  upstream  Reynolds  number  for  !:illin?T:v 
Yellow  dye  and  water  in  a  30°  branch  of  circular  cross  section  (area 
ratio  =  1.0).   Between  the  flow  ratios  of  1.0  and  0.5,  the  dye  and 
water  differ  greatly  in  Reynolds  number  but  differ  little  in  momentum 
per  unit  volume.   Note  that  the  water  also  exhibits  wall  attachment 
with  complete  main-stem  cutoff  occurring  for  Re^molds  numbers  below 
250. 
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numbers  of  dye  and  water  were  equal,  the  dye  must  also  have  had  a 
momentum  per  unit  volume  from  ten  to  fifteen  times  that  of  water.   One 
would  expect,  therefore,  that  the  greater  momentum  per  unit  volume  of 
the  entering  dye  would  carry  the  dye  preferentially  into  the  main  down- 
stream stem  resulting  in  a  lower  flow  ratio  for  the  dye  at  any  given 
Reynolds  number.   This  is  the  case  according  to  Figure  14. 

The  role  of  upstream  momentum  per  unit  volume  becomes  even 
clearer  if  one  compares  the  two  curves  in  Figure  14  at  the  same  flow 
ratio  rather  than  at  the  same  Reynolds  number. 

As  an  example,  imagine  a  horizontal  line  intersecting  the  two 
curves  at  a  flow  ratio  of  0.85.   At  this  flow  ratio,  the  corresponding 
upstream  Reynolds  numbers  are  50  and  800  for  the  dye  and  water  respec- 
tively.  That  is. 


(2X1.)       =800  (75) 

y   water 


(2£fi_.),    =  50  (76) 

^  y   '^dye 


Since  the  characteristic  length  D  is  the  same  in  each  case,  the  ratio 

of  upstream  momentum  per  unit  volume  of  water  to  that  for  the  dye  is 
given  by 


u  V^ 

((fT  ) water  ^   800   ^ water  ^  ^^     yater  ^^^^ 

(p  v)dye        50   ^'dye         ^  dye 


According  to  the  viscosity  data  given  in -Figure  6,  the  apparent  viscosity 
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of  the  dye  at  a  Reynolds  number  of  50  is  approximately  16  centipoises. 
Since  the  viscosity  of  water  may  be  taken  constant  at  1.0  centipoise, 
Equation  (77)  indicates  that  at  the  flow  ratio  of  0.85  in  Figure  lA  the 
upstream  momentum  per  unit  volume  of  the  two  liquids  is  the  same,  even 
though  their  Reynolds  numbers  differ  by  a  factor  of  16. 

From  the  above  calculation,  it  would  appear  that  the  momentum 
per  unit  volume  of  the  incoming  stream  is  the  primary  factor  determining 
flow  distribution  in  branches.   Unfortunately,  the  comparison  obtained 
for  other  flow  ratios  is  not  as  satisfactory.   For  example,  the  same 
calculation  for  a  flow  ratio  of  0.5  in  Figure  14  shows  that  the  upstream 
momentum  per  unit  volume  of  the  water  is  only  60%  of  that  for  the  Milling- 
Yellow  dye.   It  would  seem  reasonable  to  attribute  such  discrepancies, 
at  least  in  part,  to  the  presence  of  secondary  flows  at  the  higher  Reynolds 
numbers.   However,  even  with  these  discrepancies,  branching  flow  ratios 
for  different  fluids  show  a  much  stronger  correlation  with  upstream 

momentum  per  unit  volume  than  with  upstream  Reynolds  number. 

25 
Krovetz,  Crowe  and  Santi   have  sho\sm  that,  regardless  of  branching 

angle,  branches  in  general  serve  to  disturb  an  entering  flow  by  creating 

spiral  secondary  flows  in  the  branch  arms.   Since  these  disturbances  are 

also  indicated  by  a  corresponding  distortion  of  flow-birefringence 

patterns,  we  shall  characterize  the  entrance-flow  region  for  branches  by 

defining  a  bvanoh  entrance  length   for  each  arm  as  the  distance  downstream 

from  the  intersection  of  the  branch  axes  that  the  distorted  birefringence 

pattern  persists.   As  we  shall  show  in  Part  C,  it  is  possible  by  using 
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frame-by-frame  analysis  to  trace  the  propagation  of  such  disturbances 
in  time  for  nonsteady  flow  as  well. 

Figure  15  shows  the  entrance-length  data  obtained  for  a  series 

of  steady  flows  through  the  30°  branch  of  Figure  12.   Here,  the  advice 

2 
of  Atabek  has  been  followed  inasmuch  as  downstream  Reynolds  number 

has  been  used  as  the  independent  variable  rather  than  the  upstream 
Reynolds  number  used  by  Knox.^^  The  generally  accepted  correlation^'^^ 
for  entrance  flow  into  cylindrical  tubes  has  also  been  given  in  Figure 
15  for  comparison. 

It  is  apparent  that  for  Reynolds  numbers  less  than  10,  the 
entrance  lengths  in  both  arms  of  the  branch  increase  at  a  much  sharper 
rate  than  does  the  corresponding  entrance  length  in  a  straight  tube. 
This  is  undoubtedly  due  to  the  asymmetry  of  the  branch.   In  this  Reynolds- 
number  range,  the  isochromatic  patterns  show  a  smooth  transition  of 
bands  from  the  upstream  segment  into  the  arms,  and  no  band  displacement 
is  seen.   This  indicates  that  boundary- layer  separation  has  not  yet 
occurred.   Between  the  Reynolds  numbers  of  10  and  50,  both  branch 
entrance  lengths  show  no  increase,  which  indicates  that  though  more 
kinetic  energy  enters  the  branch,  this  energy  is  consumed  in  building  an 
increased  disturbance  at  the  bifurcation  rather  than  in  propagating  the 
existing  disturbance  downstream.   The  data  of  Figure  15  together  with 
the  birefringence  films  show  that  this  behavior  represents  the  onset  of 
boundary-layer  separation.   At  still  higher  Reynolds  numbers,  these 
disturbances  propagate  downstream,  but  the  fact  that  the  branch  entrance 
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Figure  15.   Branch  entrance  lengths  versus  do^mstream  Reynolds  number 
for  steady  flovr  through  a  30°  branch  of  circular  cross  section  (area 
ratio  =  1.0). 
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lengths  taper  off  shows  that  some  of  the  additional  kinetic  energy  is 
spent  in  creating  the  secondary  flows  seen  in  the  films  for  this  flow 
range.   Here  the  side-arm  entrance  lengths  are  considerably  greater 
than  those  for  the  main  stem  at  a  given  Reynolds  number,  which  agrees 
with  the  results  obtained  by  Knox.^^ 

Basically,  Figure  15  shows  that  there  are  three  phases  associated 
with  the  propagation  of  flow  disturbances  in  a  30°  branch:   1.   the 
propagation  downstream  of  an  initial  disturbance  associated  with  the 
readjustment  of  the  incoming  laminar  profile,   2.   the  onset  of  boundary- 
layer  separation  with  a  corresponding  intensification  of  disturbed  flow 
at  the  bifurcation  and  3.   the  propagation  of  this  latter  disturbance 
downstream  with  a  corresponding  development  of  secondary  flow  in  the 
arms  of  the  branch.   It  will  be  seen  shortly  that  similar  entrance-length 
patterns  are  also  found  for  flows  in  60°  and  90°  branches. 

Prior  to  investigating  branches  having  branch  angles  other  than 
30°,  it  was  felt  that  the  birefringence  patterns  and  flow-distribution 
characteristics  obtained  for  the  circular  branch  of  Figure  12  should  be 
compared  with  those  obtained  for  an  equivalent  branch  having  a  rectangular 
cross  section.   The  reason  for  this  was  twofold:   first,  to  determine 
whether  qualitative  differences  in  the  birefringence  patterns  would 
occur  owing  to  any  optical  distortion  not  accounted  for  in  the  analysis 
of  Part  A,  and,  second,   to  test  whether  a  significant  change  in  internal 
geometry  without  a  change  in  branch  angle  would  produce  corresponding 
differences  in  the  dependence  of  branch  flow  ratio  upon  upstream 
Reynolds  number. 
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Figure  16  shows  the  Isochromatic-band  distribution  obtained  with 
the  30°  rectangular  branch  (test  section  7)  for  an  upstream  Reynolds 
number  of  120,  which  is  comparable  to  the  value  of  116  for  the  30°  cir- 
cular branch  shown  in  Figure  12.   In  both  experiments,  the  incident 
light  was  circularly  polarized  and  the  area  ratio  was  unity  (the  side  arm  and 
main  downstream  stem  of  the  rectangular  branch  each  had  a  width  of  3.00 
mm  and  a  depth  of  2.00  mm). 

Qualitatively,  the  band  distributions  in  Figures  12  and  16  are 
similar,  and  the  films  indicate  that  this  similarity  in  birefringence 
patterns  prevails  over  the  entire  range  of  Reynolds  numbers  studied.  The 
bands  in  each  branch  are  displaced  toward  the  inner  walls  at  the  bifurca- 
tion, indicating  increased  shear  in  this  area,  but  neither  branch  shows 
increased  shear  at  any  other  point.   A  well-defined  stagnation  point  is 
visible  along  the  outer  side-arm  wall  in  the  rectangular  branch  (Figure 
16),  and  it  is  seen  that  the  isochroraatic  bands  tend  to  distribute  them- 
selves concavely  about  this  point.   Though  the  corresponding  point  in 

Figure  12  is  not  as  noticeable,  both  photographs  show  that  this  is  a 
» 
region  of  separated  flow.   Less  well-defined  stagnation  points  can  also 

be  seen  in  these  branches  on  the  outer  wall  of  the  main  stem  directly 

across  from  the  point  of  bifurcation. 

Figure  17  shows  the  flow-distribution  characteristics  of  the 

circular  and  rectangular  30°  branches.   It  is  apparent  from,  these  data 

that  the  marked  change  in  internal  geometry  has  little  effect  on  the 

flow  ratio  at  a  given  Reynolds  number.   The  flow-ratio  data  for  the 
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rectangular  branch  do  tend  to  fall  slightly  higher  than  those  for  the 
circular  branch,  but  this  is  probably  due  to  the  fact  that  this  branch 
possessed  a  smoother  transition  from  the  upstream  segment  to  the  side 


arm. 


Figure  18  shows  an  isochromatic-band  distribution  obtained  for 
a  low-Reynolds-number  flow  through  a  circular  60°  branch.   The  upstream 
Reynolds  number  for  this  experiment  was  10,  and  the  wall  stress  was 
75  dynes/cm2.   The  incident  light  was  circularly  polarized. 

Unlike  the  patterns  seen  at  higher  Reynolds  numbers,  this 
birefringence  pattern  is  characterized  by  a  very  smooth  and  symmetric 
transition  from  the  upstream  segment  into  the  downstream  branches. 
There  is  no  displacement  of  bands  toward  the  inner  walls,  and  the  entrance 
lengths  for  each  arm  are  less  than  two  diameters.   As  mentioned  previously, 
this  is  typical  of  the  patterns  seen  in  all  the  branches  at  very  low 
Reynolds  numbers.   The  lack  of  band  displacement  toward  the  inner  walls 
in  Figure  18  together  with  the  symmetrical  band  distribution  indicates 
that  the  flow  has  not  yet  separated.   As  the  upstream  Reynolds  number  is 
steadily  increased,  this  pattern  becomes  more  and  more  distorted  with 
the  bands  being  displaced  toward  the  inner  walls  such  as  seen  in  Figures 
12  and  16. 

Figure  19  shows  the  entrance-length  data  for  the  60"  branch. 
The  only  major  difference  between  these  data  and  those  obtained  for  the 
30°  branches  is  that  the  boundary-layer -separation  and  secondary-flow 
phases  begin  at  a  Reynolds  number  of  30  in  the  60°  branch  rather 
than  at  50  as  in  the  30°  branches.   The  entrance-length  data  in 
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Figure  19.   Branch  entrance  lengths  versus  downstream  Reynolds  number 
for  steady  flow  through  a  60°  branch  of  circular  cross  section  Carea 
ratio  =  1.0) . 
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length  data  in  Figure  19  show  the  same  three-phase  character  described 
for  the  30°  branch. 

Figure  20  shows  an  isochromatic-band  distribution  obtained  using 
plane-polarized  light  for  steady  flow  through  a  90°  branch  of  circular 
cross  section.   In  this  experiment,  the  upstream  Reynolds  number  was 
155,  the  corresponding  wall  stress  125  dynes/cm  and  the  plane  of  polari- 
zation of  the  incident  light  was  perpendicular  to  the  main  stem  of  the 
branch. 

This  birefringence  patterns  is  also  remarkably  different  from 
those  seen  for  lower  upstream  Reynolds  numbers.   A  stagnation  point  is 
visible  on  the  outer  wall  of  the  side  arm  directly  across  from  the 
bifurcation  point,  and,  more  importantly,  the  displacement  of  bands 
toward  the  inner  walls  shows  that  localized  areas  of  increased  shear 
occur  along  both  inner  walls  approximately  one  diameter  downstream  from 
the  bifurcation.   It  should  be  mentioned  that  similar  areas  were  also 
observed  at  the  highest  Reynolds  numbers  in  the  60°  branch  but  were  not 
seen  in  the  30°  branches  for  the  Reynolds-number  range  studied.   These 
observations  confirm  Fry's  speculation   regarding  sites  of  locally 
increased  shear  in  branches. 

The  shear  stress  acting  on  the  inner  downstream  walls  at  the  points 
of  maximum  band  displacement  in  Figure  20  can  be  estimated  from  Equation 
(74).   However,  in  Figure  20  the  first-order  green  band  adjacent  to  the 
inner  walls  does  not  extend  into  the  upstream  segment,  and  so  the  calcu- 
lation must  be  based  on  the  displacement  of  the  corresponding  first-order 
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green  band  appearing  in  the  upstream  segment.   That  band  has  a 
fractional  displacement  of  B(up)  =  0.4,  while  the  fractional  dis- 
placement of  the  downstream  band  is  3(down)  =  0.85  for  both  the  main 
stem  and  the  side  arm.   The  following  quantities  appearing  on  the  right 
side  of  Equation  (74)  may  then  be  obtained  from  Table  2: 


3(up)   =  0.4,   (l-g^'Cup))^  -  0.917  and  <t(up)   =  0.637 


2      % 

e(down)   =  0.85,   (1- B  (down))   =  0.527 


Since  the  upstream  wall  stress  for  this  calculation  is  125  dynes/cm^. 
Equation  (74)  gives  - 

T(avg,  down)   =   (0.917/0.527)  (0.637)  (125)   =  139  dynes/cm^. 

That  is,  for  an  upstream  Reynolds  number  of  155,  the  estimated  wall 
stress  at  the  downstream  positions  in  the  90°  branch  actually  exceeds 
the  upstream  wall  stress  by  approximately  11%. 

Figure  21  shows  the  entrance-length  data  for  the  90°  branch  of 
Figure  20.   Although  the  side  arm  of  this  branch  was  too  short  to  follow 
the  flow  disturbances  any  farther  than  four  diameters  downstream,  the 
data  for  the  main  stem  do  indicate  that  these  disturbances  propagate  in 
three  phases  just  as  for  the  30°  and  60°  branches.   The  data  for  the 
main  stem  of  the  90°  branch  show  a  more  pronounced  taper  at  higher 
Reynolds  numbers  than  do  those  for  the  other  branches,  which  indicates 
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Figure  21.   Branch  entrance  lengths  versus  downstream  Reynolds  number 
for  steady  flow  through  a  90°  branch  of  circular  cress  section  (area 
ratio  =  1.0) . 
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that  secondary  flow  in  the  90°  branch  is  also  more  pronounced  at  higher 
Reynolds  numbers.   This  secondary  flow  is  characterized  by  a  twisting 
of  the  birefringence  patterns  which  can  be  seen  in  the  corresponding 
photographs. 

Figures  22  and  23  respectively  show  the  band  patterns  for 
unseparated  and  separated  flow  through  a  +  45°  wye  branch  of  circular 
cross  section.   In  each  of  these  experiments,  the  incident  light  was 
plane  polarized  and  oriented  perpendicular  to  the  downstream  stem  appear- 
ing horizontally  in  the  photographs.   The  upstream  Reynolds  number  for 
the  unseparated  flow  of  Figure  22  was  10,  while  that  for  the  separated 
flow  of  Figure  23  was  100. 

As  was  the  case  for  the  60°  branch  shown  in  Figure  18,  the 
unseparated-flow  pattern  of  Figure  22  is  marked  by  a  smooth  transition 
of  bands  from  the  upstream  segment  into  the  downstream  arms,  with  corres- 
pondingly small  entrance  lengths.   On  the  other  hand,  the  birefringence 
pattern  of  Figure  23  exhibits  the  typical  band  displacement  and  twisting 
indicative  of  secondary  flow,  and  the  flow  disturbances  persist  a  number 
of  diameters  downstream.   The  twisting  of  the  birefringence  patterns  seen 
in  Figure  23  was  most  pronounced  in  the  90°  and  wye  branches  but  was  also 
seen  in  the  60°  branch  at  the  highest  Reynolds  numbers. 

Figure  24  shows  the  flow-distribution  characteristics  of  all  the 
branches  which  have  been  discussed.   Here,  branching  flow  ratio  has  been 
plotted  versus  the  upstream  Reynolds  number,  v/ith  branching  angle  as  a 
parameter. 
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These  data  show  clearly  that  at  any  given  Reynolds  number  the 
flow  ratio  of  single-side-arm  branches  is  only  slightly  reduced  by 
increases  in  branching  angle.   If  the  data  of  Figure  24  are  compared 
with  those  of  Figure  13,  it  is  apparent  that  variations  in  area  ratio 
produce  much  greater  variations  in  flow  ratio  than  do  variations  in 
branching  angle  for  the  same  Reynolds  number.   Figure  24  also  shows 
that  in  the  low-Reynolds-number  range,  wall  attachment  can  cause  the 
flow  ratio  to  exceed  unity  even  for  a  60°  branch.   This  effect  was  not 
observed  for  the  90°  branch.   In  contrast  to  the  single-side-arm  branches, 
the  wye  branch  showed  an  equal  distribution  of  flow  to  the  downstream 
arms  over  the  entire  Reynolds-number  range  studied. 

Finally,  the  results  of  the  above  experiments  suggest  that  there 
are  two  basic  physical  mechanisms  for  the  initiation  of  wall  trauma  and 
atherosclerotic  disease  at  sites  of  branching  in  the  arterial  system. 
According  to  the  first  mechanism, the  formation  of  atheromatous  placque 
would  be  initiated  directly   by  a  process  of  deposition  from  the  blood 
in  regions  of  local  flow  stasis.   According  to  the  second  mechanism, 
atheromatous  placque  would  be  formed  indiveatly   through  erosion  of  the 
arterial  intima  at  points  of  locally  increased  shear  followed  by  the 
fibroblastic  proliferation  and  deposition  of  mucoid  ground  substance 
described  by  Enos,  Seyer  and  Holmes.  "^^ 

The  results  of  this  study  and  that  by  Fox  and  Hugh"'-^  indicate 
that  if  atherosclerosis  results  from  the  flow-stasis  mechanism,  the  dis- 
ease would  occur  initially  along  the  outer  walls  of  branches.   On  the 
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other  hand,  if  wall  erosion  is  responsible  for  initiating  the  disease, 
early  placque  deposits  would  be  found  along  the  inner  walls. 

:^  ,  Although  several  studies  have  noted  that  atherosclerosis  shows 
a  general  predilection  for  formation  at  branches,  few  attempts  have  been 
made  to  determine  those  sites  where  the  disease  is  initiated.      This  is 
important  from  a  mechanistic  point  of  view,  since  once  deposits  are 
formed,  the  entire  geometric  character  of  the  branch  is  changed.   Of 
those  studies  reporting  on  the  distributions  of  early  lesions,  the  great 

majority  concur  with  the  findings  of  Texon^^  ^^^  Texon,  Imparato  and 

49 
Lord   who  have  shown  experimentally  that  the  point  of  bifurcation  in 

branches  is  usually  remarkably  free  of  disease,  while  the  medial-crotch 
zones  or  inner  walls  usually  show  early  involvement.   Their  results 
suggest  that  the  wall-erosion  mechanism  is  involved. 

Further  evidence  for  the  shear-stress  mechanism  is  offered  by  the 
results  of  Enos,  Beyer  and  Holmes. ^3  Reporting  a  study  of  300  coronary 
arteries  taken  from  soldiers  killed  in  Korea,  they  disclosed  that  the 
earliest  lesions  were  distributed  in  spiral  fashion  and  that  the  medial- 
crotch  walls  were  involved  a  short  distance  downstream  from  the  point  of 
bifurcation.   The  spiral  patterns  depicted  by  them  are  remarkably  similar 
to  the  twisted  birefringence  patterns  seen  at  higher  Reynolds  numbers 
in  the  present  study. 

The  question  naturally  arises  as  to  whether  shear  stresses 
occurring  in  the  cardiovascular  system  are  of  sufficient  magnitude  to 
erode  the  vessel  walls,  even  at  sites  of  locally  increased  shear.   Fryl6 
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has  determined  experimentally  that  a  shear  stress  of  379  +  85 (sd) 
dynes/cm  is  sufficient  to  erode  the  arterial  wall  in  dogs  with  exposures 
for  as  little  as  one  hour.   However,  a  similar  "critical"  value  is  not 
available  for  human  arterial  tissue. 

A  rough  estimate  of  the  maximum  wall  stress  in  any  artery  can  be 
obtained  by  neglecting  the  contribution  from  time-varying  flow  and 
assuming  that  this  stress  is  given  by  the  maximum  value  of  the  quantity 
_  (£i,)  £;  .   As  an  example,  consider  the  calculation  for  the  femoral  artery 

o2   / 

of  a  dog.  McDonald^  has  measured  the  pressure  gradient  for  this  artery, 
and  according  to  his  data, the  typical  maximum  value  of  -  (--)  is 
3  mm  Hg/cm.   The  typical  inside  radius  for  this  artery  is  0.15  cm.   Substi- 
tution of  these  quantities  above  and  converting  units  gives  an  estimated 
maximum  vrall  stress  of  300  dynes/cm^.   Since  the  stress  at  points  of 
increased  shear  would  exceed  this  value,  it  would  appear  that  at  least 
during  part  of  the  cardiac  cycle  the  wall  stress  at  these  points  would  be 
well  within  the  critical  range  reported  by  Fry.   From  the  pathologic  data 
discussed  above,  it  would  seem  reasonable  to  say  that  an  analogous  situation 
applies  to  branches  in  the  human  arterial  tree. 

C .   Pulsati le-Flov7  Experiments 

The  purpose  of  the  pulsatile-flow  experiments  was  twofold:   First, 
a  comparison  was  sought  between  the  qualitative  results  for  steady  branching 
flows  and  those  for  a  type  of  flow  more  representative  of  that  in  the 
arterial  system.   Second,  it  vjas  desired  to  investigate  the  possibility  of 
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using  flow  birefringence  to  study  quantitatively  the  propagation  of 
time-varying  disturbances. 

As  mentioned  previously,  the  pump  used  in  these  experiments  had 
the  output  flow  wave  shown  in  Figure  4,  and  experimental  data  were 
obtained  from  frame-by-frame  analysis  of  the  motion-picture  films. 

For  each  branch  studied,  no  qualitative  differences  could  be 
found  between  the  pulsatile- and  steady-flow  birefringence  patterns  for 
corresponding  upstream  flow  rates.   Typically,  as  the  initial  part  of 
the  flow  pulse  entered  the  test  sections,  the  optical  field  brightened, 
and  symmetrical  band  distributions  appeared.   These  were  similar  to  the 
distributions  shovm  in  Figures  18  and  22  for  low-Re>Tiolds'-number  steady 
flows.   Successive  frames  revealed  that  the  bands  were  then  rapidly 
displaced  toward  the  inner  walls  of  the  branches  as  regions  of  flow 
stasis  and  boundary-layer  separation  were  evidenced  along  the  outer  walls. 
As  was  the  case  for  steady  flows,  the  regions  of  flow  stasis  appeared 
directly  across  the  lumen  from  the  point  of  bifurcation.   For  those 
branches  having  branching  angles  greater  than  30°,  the  patterns  of  dis- 
placed  bands  were  followed  by  patterns  of  twisted  bands  similar  to  that 
shown  in  Figure  23.   Again,  this  behavior  was  most  pronounced  in  the  90° 
and  wye  branches.   The  succession  of  optical  phenomena  was  exactly  reversed 
for  the  descending  portion  of  the  flow  pulse. 

Significantly,  it  was  found  that  the  points  of  maximum  band 
displacement  did  not  move  with  the  incoming  flow  pulse  but  remained  cen- 
tered approximately  one  diameter  downstream  from  the  point  of  bifurcation. 
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The  regions  of  flow  stagnation  also  remained  stationary  throughout  the 
pulse  cycle.   There  was  no  evidence  to  indicate  that  any  disturbances 
were  propagated  back  upstream. 

Figures  25,  26  and  27  summarize  the  pulsatile-flow  entrance-length 
data  for  all  the  branches  which  have  been  discussed.   In  these  figures, 
branch  entrance  lengths  have  been  plotted  versus  flow-pulse  phase  angle 
(time),  and  the  flow  pulse  itself  (unlabeled  dotted  curve)  has  been 
superimposed  on  the  data  as  a  phase  reference.   In  Figure  25,  data  for 
30°  branches  of  circular  cross  section  are  compared,  with  area  ratio  as 
a  parameter.   In  Figure  26,  the  comparison  is  among  circular  branches 
having  an  area  ratio  of  unity,  with  branching  angle  as  a  parameter.   The 
data  of  Figure  27  apply  to  branches  having  an  area  ratio  of  unity,  with 
internal  geometry  and  type  of  branch  as  parameters. 

It  is  apparent  from  Figures  25  and  26  that  area  ratio  influences 
the  propagation  of  flow  disturbances  much  more  than  does  branching  angle. 
In  particular.  Figure  25  shows  that  when  the  area  ratio  was  decreased 
from  unity,  side-arm  entrance  length  decreased  and  remained  constant  over 
more  of  the  pulse  cycle,  while  that  for  the  main  stem  increased  and  tended 
to  follow  the  input  flow  wave  more  closely.   On  the  other  hand.  Figure  26 
shows  that  when  the  branching  angle  was  increased,  side-arm  entrance 
length  decreased  only  slightly  over  the  pulse  cycle,  while  that  for  the 
main  stem  increased  slightly.   These  observations  are  consistent  with  the 
fact  that  flow  ratio  decreases  markedly  with  decreasing  area  ratio  but 
is  relatively  insensitive  to  changes  in  branching  angle. 
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Figure  27  shows  that  changing  the  internal  branch  geometry  from 
a  circular  to  a  rectangular  cross  section  increased  the  entrance  lengths 
for  flow  disturbances.   This  could  be  expected  to  occur  for  any  branch 
in  which  the  cross  sections  of  the  arms  are  not  angularly  symmetric.   It 
is  also  seen  that  flow  disturbances  propagated  equal  distances  in  the 
arms  of  the  wye  branch  over  the  whole  pulse  cycle,  and  that  these  distur- 
bances propagated  farther  than  in  circular  branches  of  the  single- side-arm 
type.   This  may  be  attributed  to  the  fact  that  secondary  flow  was  more 
pronounced  in  the  wye  branch,  particularly  at  lower  flow  rates. 


SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  technique  of  flow  birefringence  has  been  used  to  investigate 
steady  and  pulsatile  flows  through  several  model  branches  simulating 
those  of  the  arterial  system.   Qualitative  and  quantitative  results 
concerning  both  the  initiation  and  propagation  of  branching-flow  distur- 
bances and  the  local  distribution  of  shearing  stress  for  such  flows  have 
been  presented.   In  addition,  flow-distribution  data  for  steady  branching 
flows  have  been  discussed. 

The  results  of  this  study  support  the  following  conclusions: 
1.   Regions  of  locally  Increased  shear,  where  stresses  can  exceed 
the  upstream  wall  stress,  occur  along  the  inner  walls  of  single- 
side-arm  branches  when  the  following  conditions  are  met:   1. 
the  upstream  Reynolds  number  exceeds  150,   2.   the  branching 
angle  is  greater  than  30°  and   3.   the  downstream  anns  of  the 
branch  are  of  approximately  the  same  diameter.   When  the  flow 
is  pulsatile,  these  regions  do  not  move  with  the  pulse  but 
remain  centered  approximately  one  diameter  downstream  from  the 
point  of  bifurcation.   These  conclusions  are  supported  by 
photographic  evidence  showing  marked  lateral  displacement  of 
isochromatic  bands  and  by  calculations  based  on  flow-birefrin- 
gence theory. 
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Regions  of  local  flow  stasis  and  boundary-layer  separation 
occur  along  the  outer  walls  of  single-side-arm  branches  for 
upstream  Rejmolds  numbers  between  20  and  50  depending  on  the 
branching  angle.   When  the  flow  is  pulsatile,  these  regions 
do  not  move  with  the  pulse  but  remain  centered  directly 
across  the  lumen  from  the  point  of  bifurcation.   These  conclu- 
sions are  supported  by  photographic  evidence  showing  stagnation 
points  characterized  by  concave  distributions  of  isochromatic 
bands. 

On  the  basis  of  conclusions  1  and  2  above  and  from  data  reported 
in  the  medical  literature, ^^'^^'22, ^7, 48, 49, 50, 54, 55  .^  .^  ^^^_ 

eluded  that  locally  increased  shear  stress  on  the  arterial 
intima  at  sites  of  branching  represents  the  most  probable  phy- 
sical mechanism  for  initiating  the  wall  trauma  associated  with 
atherosclerotic  disease.   This  contention  is  supported  by  the 
fact  that  the  maximum  wall  stress  estimated  from  McDonald's 
pressure-gradient  data^^  f^j.  ^he  femoral  artery  of  a  dog  falls 
within  the  "critical"  range  reported  by  Fry.-*"^ 
Flow  birefringence  is  a  convenient  and  useful  technique  for 
studying  the  initiation  and  propagation  of  flow  disturbances  in 
any  geometry.   Experiments  with  flows  around  needles  and  through 
branches  have  shoxm  that  flow  disturbances  can  be  characterized 
by  entrance  lengths  based  on  distorted  birefringence  patterns. 
In  particular,  it  has  been  shoxm  that  needles  mounted  flush 
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with  the  inner  wall  of  a  tube  create  no  disturbance  to  flow. 

5.  As  the  upstream  Reynolds  number  is  increased  from  zero,  flow 
disturbances  in  branches  are  generated  and  propagated  by 
three  processes:   1.   the  readjustment  of  the  initial  laminar 
velocity  profile,   2.   the  onset  of  boundary-layer  separation 
and  3.   the  onset  of  secondary  flow.   This  conclusion  is 
supported  by  entrance-length  data  for  both  steady  and  pulsatile 
flows. 

6.  For  a  given  fluid  and  geometric  configuration,  the  flow 
distribution  in  branches  is  determined  dynamically  by  upstream 
fluid  momentum  per  unit  volume  and  not  by  upstream  Reynolds 
number.   Geometrically,  the  branch  area  ratio  is  a  much  more 
significant  parameter  than  is  branching  angle.   These  conclu- 
sions are  supported  by  experiments  with  Milling-Yellow  dye  and 
water  in  the  upstream-Reynolds-number  range  from  0  to  2100. 

7.   Wall  attachment  of  flow  occurs  in  branches  of  the  single-side- 
arm  type  to  the  extent  that  side-arm  flow  will  exceed  that  in 
the  main  downstream  stem  providing  that  1.   the  upstream 
Reynolds  number  is  below  a  certain  value  which  is  dependent  on 
the  fluid,   2.   the  branching  angle  is  less  than  60°  and  3. 
the  dow-astream  arms  of  the  branch  are  of  approximately  the  same 
diameter. 

Several  recomm.endations  can  be  made  on  the  basis  of  the  above  results 
and  conclusions.   In  particular: 
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1.  The  effect  of  artificially  changing  the  pressure  pulse  in  vivo 
upon  the  incidence  and  distribution  of  atherosclerotic  lesions 
should  be  investigated  thoroughly.   Specifically,  if  the  shear- 
stress  mechanism  is  responsible  for  initiating  this  disease,  it 
would  be  expected  that  drugs  effecting  reductions  in  pressure 
gradient  while  maintaining  the  same  net  flow  would  reduce  both 
the  peak  wall  stress  and  l\\&   incidence  of  atherosclerotic  lesions 
at  branches  in  the  arterial  system. 

2.  The  use  of  micro  probes  employing  heated  thermistors,  such  as 
that  developed  by  Ling  and  Atabek,^^  should  be  investigated 
for  measuring  velocity  gradients  near  the  walls  of  branches. 
Such  measurements  could  provide  an  independent  check  on  the 
presence  of  locally  increased  shear  stress  and  could  be  made 
■in  vivo. 

3.  The  measurement  of  pressure  distributions  for  branching  flows 
would  be  a  logical  extension  of  the  present  study.   Such  experi- 
ments should  include  the  measurement  of  angular  pressure  varia- 
tions as  well  as  axial  variations.   These  results  would  be  of 
particular  importance  in  the  field  of  fluidics. 

4.  Further  studies  of  the  fluid  dynamics  of  branching  should 
incorporate  digital  computation  of  the  pressure  and  velocity 
distributiops  for  both  two-dimensional  and  three-dimensional 
branching  flows.   The  results  of  such  computations  could  not 
only  be  compared  with  the  flow-distribution  data  of  the 
present  study  but  would  provide  greater  insight  into  the 
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problem  of  branching  flows  in  general.   The  Marker  and  Cell 

53 
technique  developed  by  Welch,  Harlow,  Shannon  and  Daly   is 

ideally  suited  for  this  purpose.        '  ,. 

Both  theoretical  and  experimental  studies  of  flow  birefringence 

should  be  continued.   At  the  present  time,  laser  light  appears 

to  offer  the  most  promising  Investigative  tool  for  such 

studies. 
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APPENDIX  A 
THE  STEWART-HAMILTON  TECHNIQUE 

In  this  section  we  present  a  brief  discussion  of  the  Stewart- 

19  57 
Hamilton  indicator-dilution  technique  '    and  derive  the  commonly^  used 

working  equation  for  determining  cardiac  output  by  this  technique.   We 
demonstrate  that  this  equation  applies  only  to  a  sampling  catheter  under 
three  conditions;  namely  that  the  indicator  used  is  well  mixed  with  the 
blood  between  injection  and  sampling  sites,  that  blood  withdrawal  rate 
through  the  sampling  catheter  is  constant  and  that  recirculation  is 
eliminated.   The  first  assumption  bears  a  direct  relation  to  results  of 
the  present  study  regarding  secondary  flow  at  branches.   The  second 
assumption,  though  implied  in  the  literature,  is  seldom  stated;  and  so 
a  distinction  is  made  here  between  the  more  general  equation  including 
variable  withdrawal  and  the  special  case  used  in  practice. 

Consider  a  system  of  branching  tubes  stemming  initially  from  some 
source  trunk  such  as  indicated  in  Figure  28  below.   The  object  of  the 
Stewart-Hamilton  technique  is  to  determine  the  flow  rate  at  point  A 
(the  source)  by  rapidly  injecting  a  kno^ra  quantity  of  indicator  (say 
indocyanine-green  dye)  at  this  point  and  then  sampling  the  concentration 
of  this  indicator  as  a  function  of  time  at  an  arbitrary   point  i  down- 
stream.  It  is  shown  belov  that  the  time  integral  of  this  function  may 
be  used  to  evaluate  the  flow  rate  at  point  A.   We  wish  to  determine 
those  conditions  under  which  such  a  measurement  accurately  reflects  the 
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flow  rate  at  point  A. 
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Figure  28.   The  Stewart-Hamilton  experiment,  illustrating  a  typical 
indicator-dilution  curve.   Indicator  is  injected  at  point  A  and 
sampled  at  point  i.   Indicator  concentration  at  point  i  as  a  function 
of  time  is  shown  as  the  solid  curve  at  the  right. 

In  practice,  indicator  is  injected  and  sampled  using  cardiac 

catheters,  and  point  A,  the  injection  site,  is  located  in  or  near  the 

heart.   This  assures  that  the  resulting  flow  calculation  is  a  reflection 

of  cardiac  output  as  opposed,  for  instance,  to  flow  in  a  major  artery. 

There  are  some  complications  resulting  from  the  fact  that  the 

cardiovascular  system  is  a  closed-loop  system.   First,  since  several 

injections  may  be  made  during  the  course  of  a  single  cardiac-catheteri- 

zation  procedure,  account  must  be  made  for  the  presence  of  any  indicator 
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due  to  previous  injections.   This  indicator  is  reflected  in  the  base-line 
concentration  level  Cq  shown,  in  the  figure.   The  situation  is  further 
complicated  by  the  fact  that  conimonly  used  indicators  are  degraded  by 
certain  organs  of  the  body  and  consequently  exhibit  a  half -life.   This 
half-life,  approximately  twenty  minutes,  is  considerably  greater  than 
the  sampling  time  tg,  which  is  of  the  order  of  twenty  seconds;  so  the 
associated  base-line  shifting  does  not  lead  to  appreciable  error  in 
cardiac-output  calculations.   Secondly,  typical  indicator-dilution  curves, 
such  as  that  pictured  in  the  figure,  exhibit  a  recirculation  hump  owing 
to  indicator  having  traversed  the  loop  and  making  a  second  pass.   Thus, 
for  accurate  results, these  curves  must  be  extrapolated  in  such  a  way  as 
to  disallow  any  contribution  from  recirculation.   An  exponential  extra- 
polation is  a  common  and  accurate  way  of  doing  this. 

Fortunately,  these  complications  can  be  easily  accounted  for 
during  an  actual  procedure;  and,  in  fact,  an  analog-computer  program 
has  been  developed   which  automatically  does  this  as  well  as  providing 
immediate  cardiac  outputs. 

Referring  again  to  Figure  28,  we  now  wish  to  derive  the  working 
equation  of  the  Stewart-Hamilton  technique  and  to  consider  the  assumptions 
involved  in  obtaining  this  relation. 

It  is  assumed  in  this  derivation  that  indicator  concentration  at 
point  i  is  measured  along  the  extrapolated  dilution  curve.   This  is  tanta- 
mount to  assuming  that  injected  indicator  passing  point  i  makes  only  a 
single  pass,  or  that  recirculation  has  been  eliminated. 
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Let  us  now  assume  that  a  known  amount  of  indicator  (say  A  mg) 
is  rapidly  injected  at  point  A.   This  indicator  simultaneously  mixes 
with  the  blood  flov/ing  past  point  A,  diffuses  and  is  convected  down- 
stream.  Since  axial  diffusion  is  negligible  compared  to  axial  convec- 
tion, the  amount  of  indicator  which  passes  any  arbitrary  downstream 
point  i  in  the  time  interval  dt  may  be  expressed  as 

dA^  =  QiCidt  (Al) 

where:         Q-  =  volumetric  flow  rate  at  point  i,  (1/sec) 

Cj   =  indicator  concentration  at  point  i,  (mg/1) 

If  time  is  measured  starting  with  the  injection  at  time  zero,  then 

the  total  amount  of  indicator  which  passes  point  i  from  t  =  0  to  t  =  t„  is 

ft, 
A,  =  1  0,C,dt  (A2) 


.,  =  l^^c,. 


where  the  sampling  time  t  is  sufficiently  great  to  allow  all  indicator 
pursuing  this  particular  route  to  pass  point  i.   It  should  be  mentioned 

that  an  equation  analogous  to  Equation  (A2)  has  been  given  previously  by 

12 
Cropp  and  Burton.    However,  their  expression  replaces  the  left  side  of 

Equation  (A2)  with  the  total  amount  of  indicator  injected  and  is  therefore 

only  applicable  to  a  non-branching  system.   They  did  not  discuss  a 

branched  system,  which  is  clearly  of  much  more  interest. 

We  now  make  the  first  major  assumption.   It  is  assumed  that  the 

indicator  injected  at  point  A  becomes  well  mixed  over  the  imaginary 
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surface  perpendicular  to  the  flow  streamlines  at  this  point  and 
subsequently  remains  well  mixed  over  all  such  surfaces  throughout  the 
flow  field  between  injection  and  sampling  sites.   This  is  to  say  that 
over  any  surface  in  the  flow  field  which  is  everywhere  perpendicular 
to  the  flow  streamlines,  indicator  concentration  depends  only  on  time 
and  not  on  position.   In  particular,  this  assumption  is  meant  to  apply 
at  branches.   For  the  special  case  of  parallel  flow  in  a  circular  tube, 
this  is  equivalent  to  neglecting  any  angular  or  radial  concentration 
gradients  but  not  axial  ones. 

The  above  assumption  might  be  justified  in  practice  due  to  a 
combination  of  factors  such  as  mixing  produced  in  the  ventricles  of  the 
heart,  secondary  flow  at  branches  and  a  resultant  "flat"  velocity  profile 
due  to  the  pulsatile  nature  of  the  flow.   At  best,  however,  it  represents 
an  approximation.   A  study  by  Bassingthwaighte,  Edwards  and  Wood^  has 
shown  that  there  is  good  agreement  between  cardiac  outputs  determined 
simultaneously  from  central  and  peripheral  sites  in  dogs.   Their  results 
lend  support  to  the  assumption  of  well-mixed  indicator,  but  only  for  the 
larger  vessels. 

Note  that  axial  concentration  gradients  still  exist  in  the 
cardiovascular  system,  as  is  evident  from  the  nature  of  a  typical  dilution 
curve;  thus,  the  assumption  of  well-mixed  indicator  is  meant  only   in  the 
sense  discussed  above.   This  point  has  not  been  clearly  made  in  the 
literature. 

Since  the  assumption  of  well-mixed  indicator  infers  that  indicator 
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Injected  at  point  A  uniformly  tags  the  blood  over  the  above-mentioned 
Imaginary  surfaces  from  this  point  on,  it  is  Implied  that  the  total 
amount  of  indicator  entering  any  downstream  branch  from  t  =  0  to  t  =  t 
is  directly  proportional  to  the  fraction  of  the  total  flow  past  A  which 
enters  that  branch  in  the  same  interval.   Stating  this  mathenatically, 
the  assumption  of  well-mixed  indicator  implies  the  following: 

ts 


f  Qi^^ 

-2- 
t. 


Ag  =  K  — _  (A3) 


s 

dt 


0 


where:   A^  =  total  amount  of  indicator  (mg)  which  enters  the 
downstream  branch  in  which  point  i  is  located 
in  the  time  interval  from  t  =  0  to  t  =  tg 

K   =  constant  of  proportionality 

Qi  =  volumetric  flow  rate  at  point  i,  (  1/sec  ) 

Q^  =  volumetric  flow  rate  at  point  A,  (  1/sec  ) 

To  evaluate  the  constant  of  proportionality  K,  let  the  experiment 
be  conducted  in  an  unbranched  tube.    For  this  case,  Qj  =  Q^  and  Ag  = 
A,  the  total  amount  of  indicator  injected.   Thus,  the  proportionality 
constant  is  equal  to  A,  and  the  total  amount  of  indicator  entering  the 
arbitrary  downstream  branch  becomes 


Ae  =  A_y (A4) 

■  s 
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Since  there  is  no  net  accumulation  of  indicator  in  this 
arbitrary  branch,  one  may  set  the  total  amount  of  indicator  which 
enters  according  to  Equation  (A4)  equal  to  the  total  amount  which 
passes  point  i  according  to  Equation  (A2)  to  obtain 


*.s  t 

s 


f  Q.dt 
•'0 


QA<it 


QiC.dt  (A5) 


If  one  defines  an  integral  mean  of  the  quantity  Q  as 

t 

Q(t)  dt 


I 


Q(mean)   =  ■ ■  (A6) 


I 


dt 
0 


then  the  left  side  of  Equation  (A5)  can  be  rewritten  in  terms  of  mean 
flow  rates  for  the  interval  t  =  0  to  t  =  tg  by  dividing  numerator  and 
denominator  by  tg.   Thus, 

t 
s 

Q . (mean) 

A     ^  .   I  Q.C.dt  (A7) 


I  Q.C.c 


Since  the  quantity  of  interest  in  the  Stewart-Hamilton  experiment 
is  Q.  (mean),  we  solve  Equation  (A7)  to  obtain 

A  Q£(mean) 

Q.  (mean)   = ~  (A8) 

A  t 

s. 
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This  is  a  general  form  of  the  Stewart-Hamilton  working  equation 
for  determining  cardiac  output,  and  it  applies  for  time-varying  flow 
and  concentration  at  any   downstream  sampling  site.  The  assumptions 
involved  in  this  expression  are  well-mixed  indicator  in  the  sense 
discussed  previously  and  the  absence  of  recirculation. 

It  is  evident  that  to  apply  Equation  (A8)  rigorously  to  a 
downstream  site  (such  as  the  femoral  artery)  one  would  need  to  know  not 
only  concentration  C.  but  also  flow  rate  Q.  as  functions  of  time  at  this 
site.   Unfortunately,  this  information  is  generally  unavailable  from 
catheterization  data;  and,  in  fact,  since  concentration  is  always  meas- 
ured externally  by  withdrawing  blood  through  a  sampling  catheter  and 
transducer,  neither  concentration  nor  flow  rate  is  actually  measured  at 
point  i. 

This  apparent  difficulty  may  be  overcome  in  the  following  way: 
Since  Equation  (A8)  applies  to  any  downstream  sampling  site,  it  may  be 
taken  to  apply  in  particular  at  the  transducer  site  in  the  withdrawal 
catheter  where  concentration  actually  is  measured.   If  Equation  (A8)  is 
applied  in  this  way,  then  flow  rate  Q^  becomes  the  blood  withdrawal  rate 
through  the  sampling  catheter,  which  may  be  controlled  externally,  and 
C-i  is  the  properly  scaled  output  of  the  transducer.   In  particular,  if 
blood  is  withdrawn  at  a  constant   rate,  then  Q^^  is  constant  and  identically 
equal  to  Q.(mean).   The  resulting  simplification  of  Equation  (AS)  gives 
the  working  equation  for  determining  cardiac  output  as  it  is  usually 

employed.   That  is 


^v^me^if?:^ 
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Q.(mean)   =  — (A9) 


|c,dt 


0 

It  is  important  to  note  that  Equation  (A9)  applies  to  the 
withdrawal  catheter  and  is  only  valid  for  constant  withdrawal  rate, 
in  addition  to  being  subject  to  the  assumptions  of  well-mixed  indicator 
and  no  recirculation.   If  withdrawal  rate  is  not  constant  or  if  flow 
and  concentration  measurements  actually  are  made  at  a  point  i  in  the 
cardiovascular  system,  then  the  more  general  equation  (A8)  must  be 
used. 
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Table  3 


APPARENT-VISCOSITY  AND  REYNOLDS-NUMBER  DATA 
FOR  1.4%  MILLING- YELLOW  DYE 
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